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ABSTRACT 


This  report  discusses  advances  in  the  development  of  a  multi-frequency 
modulation  (MFM)  packet  communications  system  on  an  industry  standard 
computer.  Transmitter  and  receiver  programs  are  described  that  control  vector 
signal  processors  and  data  acquisition  devices.  Further,  these  programs  encode, 
modulate,  demodulate,  and  decode  the  MFM  signal.  The  throughput  data  rate 
was  doubled,  the  encoding/decoding  process  was  implemented  in  near  real-time, 
and  a  personal  computer  plug-in  board  was  designed  and  built  to  provide  syn¬ 
chronization  between  the  transmitter  and  receiver.  This  MFM  implementation 
provided  acceptable  bit  error  rates  at  input  signal-to-noise  ratios  of  15  dB  and 
above. 
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I.  INTRODUCTION 


Multi-frequency  modulation  (MFM)  provides  many  advantages  for  typical 
digital  communication  links.  The  entire  system  is  implemented  on  industry 
standard  personal  computers  (PC).  Analog  devices,  such  as  voltage  multipliers 
and  phase-locked-loops,  are  not  required  in  the  modulation  of  the  baseband  sig¬ 
nal  or  in  the  demodulation  of  the  received  bandpass  signal.  The  frequency  range 
of  the  signal  is  limited  only  by  the  direct  memory  access  (DMA)  rate  of  the  PC 
and,  within  this  range,  is  easily  changed  under  software  control.  Also, 
guardbands  are  essentially  eliminated,  allowing  much  more  efficient  use  of  the 
available  frequency  spectrum.  Finally,  the  MFM  technique  as  implemented  here 
exclusively  uses  differential  encoding  schemes,  eliminating  the  need  for  coherent 
detection. 

The  scope  of  this  thesis  covers  the  improvement  of  an  existing  MFM  com¬ 
munication  system  developed  by  LT  Terry  K.  Gantenbein,  USN  [Ref.  1],  New 

concepts  implemented  during  this  project  include  the  following: 

•  the  use  of  vector  signal  processors  to  make  the  signal  encoding  and  decoding 
process  near  real-time, 

•  the  design  of  software  encoding  algorithms  to  double  the  data  rate  of  the 
existing  system,  and 

•  the  design  and  implementation  of  a  PC  plug-in  board  used  to  synchronize 
the  receiver  to  the  transmitter. 

The  theory  of  MFM  is  briefly  reviewed  in  Chapter  II,  the  system  design  is 
discussed  in  Chapter  III,  the  implementation  is  delineated  in  Chapter  IV,  and  the 
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analysis  of  the  system  is  presented  in  Chapter  V.  The  reader  interested  in  a  more 
thorough  discussion  of  the  properties  of  MFM  is  referred  to  Refs.  1  and  2. 
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II.  BACKGROUND  OF  MULTI-FREQUENCY  MODULATION 
A.  THE  MFM  SIGNAL  PACKET 

The  MFM  signal  packet  has  time  and  frequency  characteristics  which  give  it 
unique  properties.  The  following  terms  are  used  in  the  description  of  the  MFM 
signal  packet: 

T :  Packet  length  in  seconds 
AT :  Baud  length  in  seconds 
kx  :  Baud  length  in  number  of  samples 
4>lk  :  Phase  of  the  k'h  tone  in  the  Ith  baud 
Alk  :  Amplitude  of  the  kth  tone  in  the  Ith  baud 
L  :  Number  of  bauds  per  packet 
At :  Time  between  samples  in  seconds 
fx  =  1/A? :  Sampling  frequency  or  A/D  and  D/A  conversion  rate  in  Hz 

A f  =  \jAT :  Frequency  spacing  between  tones 
K :  Number  of  MFM  tones  in  a  baud 

These  MFM  parameters  are  depicted  in  Figure  1. 

The  MFM  signal  packet  is  composed  of  a  series  of  "bauds"  of  length  AT  . 
By  choosing  AT  =  /c,A?,  each  baud  consists  of  kz  time  samples  spaced  At  seconds 
apart.  The  vertically  stacked  bins  in  Figure  1  represent  the  frequencies  (tones) 
in  each  baud.  Each  baud  has  the  same  tones,  but  the  tones  are  encoded  with 
different  information  from  baud  to  baud.  One  property  of  MFM  is  that  the 


3 


Figure  1.  Time/frequency  parameters  of  MFN1  signal  packet 

frequencies  of  each  of  the  tones  in  a  baud  are  integer  multiples  of  A/,  thus  the 
tones  are  harmonics  with  frequency  spacing  A/  Hz  [Ref.  2:  pp.6-7].  In 
Figure  1,  tones  f2  and  y,  are  ( k2  -  A.-,)A/  Hz  apart.  The  sampling  theorem  requires 
that  the  highest  frequency  in  the  signal  must  not  exceed  one-half  the  sampling 
rate,/./2.  Thus,  in  an  MFM  signal  packet  the  highest  representable  frequency  is 
fJ2  -  A if. 

In  Figure  1,  symbol  'Ik'  is  a  single  tone  in  the  baud  with  frequency  kAf, 
amplitude  Alk,  and  phase  <f),k.  The  amplitude  and  phase  of  the  tone  are  assigned 
in  accordance  with  any  suitable  encoding  scheme.  Selection  of  a  scheme  to  rep¬ 
resent  four  bits  per  tone  is  a  focal  point  of  this  thesis  and  is  presented  in  Chapter 
III.  The  frequency  range  of  the  signal  is  assigned  by  choosing  the-harmonic  val¬ 
ues  and  k2  to  coincide  with  the  desired  passband.  Herein  lies  one  of  the 


flexibilities  of  MFM.  Given  a  sampling  frequency,  the  desired  passband  is 
selected  by  simply  assigning  the  correct  values  for  fc,  and  k2. 

The  mathematical  description  of  the  analog  MFM  signal  in  any  baud  is  given 
by  [Ref.  2:  pp.6-7] 


*,/  2—1 

X({t)  =  ^  Alk  cos(2nkAft  +  (/  -  1)AT  ^  t  <,  IAT  (1) 

k=  1 

and  the  sampled  discrete  time  representation  is 


kxl  2-1 

Xj{n)  =  ^  Alk  cos(  +  4>lk ),  0  £  n  <,  kx  -  1 .  (2) 

A=1 

B.  MFM  SIGNAL  PACKET  CONSTRUCTION  AND  MODULATION 

The  Inverse  Discrete  Fourier  Transform  (IDFT)  with  its  symmetry  properties 
is  ideally  'uited  for  the  generation  of  the  MFM  signal  packet.  A  symmetry 
property  of  the  Discrete  Fourier  Transform  (DFT)  is  the  following:  the  DFT  of 
a  /c,-point  real-valued  sequence  is  a  kx- point  complex  sequence  with  conjugate 
image  symmetry.  That  is,  the  values  in  the  second  half  of  the  complex  sequence 
will  be  the  complex  conjugates  of  their  image  values  in  the  first  half  of  the  com¬ 
plex  sequence.  This  relationship  is  diagrammed  in  Figure  2  for  kx=  16  points. 
Conversely,  the  IDFT  of  a  conjugate  symmetric  array  of  frequency  domain  val¬ 
ues  produces  an  entirely  real  sequence  of  time  domain  values.  This  property  is 
exploited  as  follows.  After  all  the  tones  have  been  encoded,  the  amplitudes  and 
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Figure  2.  IDFT  data  structure 


phases  are  converted  from  polar  to  rectangular  coordinate  form.  Next,  these 
rectangular  frequency  plane  coordinates  are  loaded  into  the  First  half  of  a  kx- point 
complex  array  and  the  complex  conjugates  of  the  values  in  the  first  half  of  the 
array  are  placed  in  the  second  half  of  the  array  at  the  image  frequencies.  The 
IDFT  is  performed  on  this  array,  yielding  kx  real  time-domain  values.  This  proc¬ 
ess  is  repeated  for  each  baud,  giving  kxL  real  samples.  At  this  point,  the  packet 
is  ready  to  be  modulated. 

Modulating  the  packet  is  simpler  than  creating  it.  Since  all  the  frequency  in¬ 
formation  was  assigned  in  the  selection  of  the  passband,  all  that  must  be  done  to 
modulate  the  packet  is  to  pass  it  through  a  digital-to-analog  (D/A)  converter  at 
clock  rate  fx.  The  block  diagram  of  the  MFM  transmitter  is  given  in  Figure  3. 
The  resultant  D/A  output  contains  the  desired  frequencies  in  the  range  initially 
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Figure  3.  Block  diagram  of  the  MFM  transmitter 


selected,  k{Af to  k2Af.  The  passband  used  in  this  thesis  models  an  acoustic  chan¬ 
nel  in  the  range  16-20  KHz;  the  approximate  MFM  parameters  are  shown  in 
Table  1  below. 


Table  1.  PARAMETERS  FOR  A  1/I5TH  SECOND  MFM  SIGNAL  PACKET 
IN  A  16-20KHZ  PASSBAND 


Baud  Iength(sec) 

AT 

1/240 

1/120 

1/60 

1/30 

1/15 

No.  of 
bauds 

L 

16 

8 

4 

2 

1 

Tone  spacing 

A/ 

240 

120 

60 

30 

15 

Lowest  harmonic 

k\ 

68 

135 

269 

537 

1073 

Lowest  tone  freq. 

fi 

16320 

16200 

16140 

16110 

16095 

Highest  harmonic 

k2 

83 

166 

332 

664 

1328 

Highest  tone  freq. 

h 

19920 

19920 

19920 

19920 

19920 

Samples  per 
baud 

kx 

256 

512 

1024 

2048 

4096 

Sampling  freq 

fx 

61440 

61440 

61440 

61440 

61440 
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Note  that  any  reference  to  the  "size"  of  a  baud  describes  both  the  number  of 
samples  in  the  baud  and  the  number  of  points  in  the  DFT  used  to  generate  the 
samples  in  the  baud.  The  reader  interested  in  voice-band  channel  applications  is 
referred  to  Salsman  [Ref.  3]. 

C.  THE  MFM  RECEIVER 

The  MFM  reception  process  is  the  reverse  of  the  MFM  transmission  process. 
After  the  receiver  is  synchronized  to  the  transmitter  (see  Chapter  III),  the  data 
is  sampled  through  an  analog-to-digital  (A/D)  converter  at  the  sampling  fre¬ 
quency  fx.  A  ^,-point  DFT  is  performed  on  kx  samples.  The  second  half  of  the 
array  of  frequency  domain  values  is  discarded  as  it  is  redundant  information. 
Next,  the  tones  are  decoded  and  the  data  is  presented  to  the  information  sink. 
This  DFT/decoding  loop  is  repeated  L  times  until  all  the  Lkx  samples  have  been 
processed.  The  MFM  receiver  block  diagram  is  shown  in  Figure  4. 


Figure  4.  MFM  receiver  block  diagram 
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III.  SYSTEM  DESIGN 


This  chapter  presents  the  modifications  made  to  the  MFM  communication 
system  developed  by  Gantenbein  [Ref.  1].  First,  the  encoding  process  will  be 
discussed  with  emphasis  given  to  signal  constellations,  differential  encoding,  and 
bit  rates.  Next,  receiver  synchronization  and  its  relationship  to  system  phase  re¬ 
sponse  is  presented.  This  section  will  also  cover  the  development  of  a  synchronizer 
PC  plug-in  board  for  the  receiver  terminal.  Finally,  the  use  of  vector  signal 
processors  in  the  computation  of  the  Fourier  transforms  will  be  reviewed. 

A.  ENCODING  AND  DECODING 

The  assignment  of  amplitudes  and  phases  to  the  MFM  tones  was  mentioned 
in  Chapter  II.  This  process  is  known  as  encoding  and  may  be  based  on  any 
realizable  signal  constellation.  Further,  in  the  MFM  research  conducted  at  NPS, 
the  data  encoded  has  been  files  of  symbols  from  the  American  Standard  Code  for 
Information  Interchange  (ASCII)  character  set.  Though  ASCII  characters  were 
used  here,  any  binary  information  source  may  be  used. 

Digital  communications  techniques  generally  fall  into  two  categories,  coherent 
and  noncoherent.  Coherent  systems  must  extract  the  carrier  phase  information 
from  the  incoming  signal  in  order  to  decode  a  received  symbol,  which  is  difficult 
and  expensive.  Noncoherent  differential  systems  do  not  attempt  to  recover  any 
carrier  phase  or  amplitude  information.  Instead,  decoding  is  performed  by  com¬ 
paring  phase  and/or  amplitude  differences  between  adjacent  received  symbols. 
Theoretically,  MFM  can  be  implemented  using  either  coherent  or  differential 
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noncoherent  demodulation.  For  the  coherent  case,  since  the  tones  are  all 
harmonically  related,  a  pilot  tone  could  be  placed  in  a  frequency  bin  in  each  baud 
several  bins  away  from  the  rest  of  the  message.  The  phase  of  this  tone  could  then 
be  recovered  using  traditional  carrier  phase  recovery  loops.  However,  this  still 
may  be  inadequate  to  determine  the  phase  of  all  the  tones  in  an  MFM  signal  in 
channels  where  there  is  significant  uncompensated  group  delay.  In  the  nonco¬ 
herent  case,  a  differential  encoding  scheme  is  used  to  encode  the  tones;  no  carrier 
phase  information  is  required. 

Though  the  differential  schemes  do  not  require  carrier  phase  recovery,  de¬ 
tection  performance  will  be  two  to  three  dB  inferior  to  a  coherent  scheme  using 
the  same  signal  constellation.  For  example,  to  achieve  a  probability  of  bit  error 
of  10~-\  the  input  signal-to-noise  ratio  ( SNRm )  for  differentially  encoded  phase- 
shift-keying  (DPSK)  is  8.5  db,  while  for  coherent  phase-shift-keying  it  is  6.5  db 
[Ref.  4:  p.  160].  This  thesis  research  studied  noncoherent  encoding/decoding 
schemes  exclusively. 

The  data  rate  (without  error  correction  coding)  in  an  MFM  system  is  deter¬ 
mined  by  the  number  of  bits  encoded  in  each  tone.  Gantenbein  [Ref.  1]  used  a 
differential  quadrature-phase-shift-keyed  (DQPSK)  modulation  constellation 
which  encoded  two  bits  into  each  tone.  The  bandwidth  efficiency,  Bv  ,  in  any 
MFM  packet  communications  system  is  given  by 


B. 


_  /  No.  bits  \  /  tones  \  ,  bauds 
'  tone  ’ '  J 


baud  second  bandwidth 


)( 


1 


)• 


(3) 


10 


For  example,  using  the  values  for  any  of  the  baud  sizes  in  Table  1,  the  band¬ 
width  efficiency  is  two  bits/sec/Hz  for  DQPSK  which  encodes  two  bits  into  each 
tone.  It  is  apparent  then  that  even  though  the  MFM  parameters  vary,  the 
bandwidth  efficiency  in  bits/sec/Hz  is  the  number  of  bits  encoded  in  each  tone. 

The  challenge  is  to  develop  a  signal  constellation  and  differential 
encoding/decoding  scheme  that  will  increase  the  bit  rate  while  giving  acceptable 
SNR  and  corresponding  probability  of  bit  error  performance.  To  this  end,  two 
separate  differential  16  quadrature-amplitude-modulation  (D16QAM) 
constellations  are  described  below. 

1.  D16QAM  constellation  one 

The  first  D16QAM  signal  constellation  studied  is  shown  in  Figure  5.  In 
this  constellation,  the  complex  plane  is  divided  into  45°  sectors  and  two  concen¬ 
tric  "magnitude''  circles.  The  inner  tone  magnitude  circle  is  labeled  as  IX  and  the 
outer  tone  magnitude  circle  (which  is  twice  the  magnitude  of  1A^  is  labeled  as 
2X.  The  resultant  complex  plane  has  16  distinct  regions,  each  of  which  is  assigned 
a  four  bit  value.  The  sectors  are  numbered  in  a  counterclockwise  direction  from 
the  positive  real  axis  starting  with  the  inner  magnitude  circle. 

Using  this  signal  constellation,  the  differential  encoding  for  one  symbol  is 
show-n  in  Figure  6.  Assume  that  the  previously  encoded  tone  in  a  baud  is  at 
position  A  with  phase  67.5°  as  seen  in  step(l).  An  ASCII  character  is  retrieved 
from  the  message  to  be  transmitted  and  the  first  four  bits  are  retained.  The  last 
four  bits  are  saved  for  the  next  encoding.  The  least  significant  (rightmost)  bit  of 
the  four  bit  symbol  is  stripped  off  and  saved,  and  the  value  of  the  remaining  three 
bits  is  multiplied  by  45°.  The  resulting  value  (in  the  example,  90°)  represents  the 
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Figure  5.  16QAM  signal  constellation 

amount  of  phase  rotation  to  add  to  the  previously  encoded  tone's  phase  to  arrive 
at  the  new  tone's  phase.  As  shown  here  in  step(3),  the  new  tone  is  assigned  a 
phase  value  of  157.5°  . 

Next,  the  magnitude  of  the  newly  encoded  tone  is  determined  using  the 
bit  that  was  originally  removed  from  the  four  bit  symbol  in  step(2).  The  magni¬ 
tude  assignment  is  much  simpler  than  the  phase  assignment.  If  the  least  signif¬ 
icant  bit  is  a  zero,  the  magnitude  of  the  new  tone  is  made  the  same  as  the 
previous  tone's  magnitude.  If,  however,  the  least  significant  bit  is  one,  the  new 
tone  is  assigned  a  magnitude  different  than  the  previous  tone's.  In  step  (2),  the 
least  significant  bit  value  is  zero  and  the  previous  tone  magnitude  was  2X,  so  the 
new  tone  magnitude  will  also  be  2A'. 
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Figure  6.  Encoding  example  with  D16QAM 

This  constellation  is  referred  to  as  the  'magnitude  encoding  scheme'  since 
the  change  of  phase  assigned  to  the  tone  is  based  exclusively  on  the  magnitude  of 
the  three  bits  from  each  symbol. 
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In  the  receiver,  decoding  is  performed  to  extract  the  tone  phase  and 
magnitude  information.  First,  two  adjacent  tones  are  multiplied  together  using 
complex  arithmetic.  The  phase  difference  A<f>  is  found  from  [Ref.  1:  p.l  1] 

gf**  -  !)*  =  gttig-j*!- 1  _  gi(<l>i-4>i- 1)  ^ 

In  the  absence  of  noise,  A <f>  will  be  an  integer  multiple  of  45°.  This  difference  is 
then  realigned  to  the  original  signal  constellation  by  rotating  it  counterclockwise 
22.5°,  and  the  bit  pattern  in  that  sector  is  chosen  as  the  First  three  bits  in  the 
symbol. 

The  last  bit  in  the  symbol  is  recovered  by  comparing  the  absolute  magni¬ 
tudes  of  the  two  adjacent  tones.  Since  the  transmitted  tones  may  have  only  one 
of  two  amplitudes,  the  decision  region  is  placed  at  1.5X ,  or  at  one-half  the  dis¬ 
tance  between  \X  and  2X.  So  if  a  tone's  magnitude  is  within  one  and  one-half  the 
magnitude  of  the  previous  tone  it  is  decoded  as  the  same  magnitude  and  the  last 
bit  is  decided  to  be  a  zero.  If  the  tone's  magnitude  is  not  within  this  range,  it  is 
decoded  as  having  a  different  magnitude  and  the  bit  is  decided  to  be  a  one. 

2.  Gray-encoded  DI6QAM  constellations 

Gray-encoded  signal  constellations  reduce  the  probability  of  bit  error  in 
the  presence  of  Additive  White  Gaussian  Noise  (AWGN).  In  a  non-Gray- 
encoded  constellation,  the  bit  patterns  in  adjacent  sectors  may  differ  in  as  many 
as  three  positions,  as  shown  in  sectors  one  and  eight  in  Figure  5. 

The  D16QAM  constellation  shown  in  Figure  7  differs  from  the  magni¬ 
tude  constellation  in  that  each  sector  differs  in  only  one  bit  position  with  respect 
to  any  adjacent  sector.  To  encode  a  tone,  the  phase  rotation  to  be  added  to  the 
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Figure  7.  Gray-encoded  16QAM  signal  constellation 

previous  tone's  phase  is  obtained  from  a  lookup  table;  there  is  no  simple  arith¬ 
metic  relationship  between  the  magnitude  represented  by  the  three  bit  pattern 
and  the  amount  of  rotation.  The  magnitude  bit  is  encoded/decoded  as  before. 

B.  SYNCHRONIZATION 

Digital  communication  links  require  different  levels  of  synchronization,  based 
on  the  particular  application  and  modulation  type.  MFM  communication  sys¬ 
tems  require  packet  synchronization  at  the  receiver.  The  receiver  PC  must  have 
some  method  for  initiating  sampling  at  its  end  of  the  channel.  This  is  accom¬ 
plished  in  the  following  manner. 

The  transmitted  packet  contains  a  synchronization  baud  (syncbaud)  ap¬ 
pended  to  its  front  end  as  shown  in  Figure  8.  The  syncbaud  is  a  baud  of  size  256, 
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Figure  8.  Synchronizer  functional  block  diagram 


generated  in  the  same  manner  as  the  other  bauds  in  the  packet.  The  receiver  has 
a  synchronization  device  that  has  as  its  primary  feature  a  128-bit  correlator.  The 
last  128  bits  of  the  syncbaud  sent  by  the  transmitter  are  also  stored  in  the  receiver 
and  used  as  the  reference  sequence  for  the  correlator.  The  correlator  is  clocked 
at  the  sampling  frequency  fx,  and  compares  the  polarities  of  the  incoming  signal 
to  the  reference  polarities  used  by  the  transmitter.  In  this  manner,  the 
synchronizer  yields  a  correlation  peak  at  the  last  sample  of  the  syncbaud.  This 
peak  triggers  a  latch  which  enables  the  sampling  clock  at  the  A/D  converter  in 
the  receiver  PC.  Thus,  sampling  commences  at  the  first  sample  of  the  first  infor¬ 
mation  bearing  baud  of  the  packet.  Gantenbein  [Ref.  1:  pp.  24,38]  implemented 
synchronization  with  a  polarity  coincidence  correlator  on  a  breadboard  circuit 
external  to  the  receiver  PC.  His  prototype  provided  adequate  performance  but 
was  too  large  to  implement  on  a  PC  plug-in  card. 

Recent  large  scale  integration  (LSI)  chip  technology  enabled  this  author  to 
perform  the  same  128  point  correlation  on  a  plug-in  card,  with  the  added  feature 
of  being  able  to  load  the  correlator  reference  sequence  from  the  receiver  PC.  An 
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overview  of  this  device  is  presented  here;  a  detailed  description  and  schematic  are 
given  in  Appendix  A. 

The  PC  card-based  synchronizer  was  developed  using  the  TRW  LSI  Pro¬ 
ducts,  Inc.,  TMC222I  correlator  chip  as  its  centerpiece.  The  TMC2221  is  a 
128-point  programmable  digital  correlator  implemented  in  a  28  pin  dual  inline 
package  (DIP). 

Two  separate  functions  occur  on  the  synchronizer  PC  board  whose  block  di¬ 
agram  is  shown  in  Figure  9.  First,  the  reference  values  are  loaded  into  the  cor¬ 
relation  reference  register  on  board  the  TMC2221  chip.  Direct  memory  access 
(DMA)  in  the  single  byte  transfer  mode  is  used  to  send  the  reference  values  from 
the  PC  memory  to  the  chip.  The  single  byte  transfer  mode  implies  that  a  device 
requesting  DMA  transfer  must  provide  a  separate  DMA  request  for  each  byte  it 
wants  placed  on  the  PC  data  bus.  To  obtain  the  128  reference  values,  the 
synchronizer  board  must  make  128  byte  transfer  requests.  The  data  transfer  ac¬ 
knowledge  signals  returned  on  the  PC  bus  by  the  DMA  controller  are  then  used 
to  control  the  reference  register  loading  in  the  TMC2221  chip. 

DMA  is  typically  used  to  provide  high  speed  data  transfers  from  or  to  PC 
memory,  usually  on  the  order  of  thousands  of  bytes  per  second  (as  is  done  in  the 
MFM  transmitter  D/A  process  to  clock  out  up  to  61440  byte  samples  per  second 
[Ref.  1,5]).  Its  use  here  was  selected  with  the  vision  that  it  would  minimize  the 
overall  circuitry  required  to  implement  the  transmitter  and  receiver  hardware  on 
a  single  PC  plug-in  card  in  the  future.  Further,  because  the  reference  values  are 
easily  changed,  different  acquisition  delays  could  be  tested,  thus  enabling 
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Figure  9.  Block  diagram  of  PC  synchronizer  board 

optimization  of  the  system  phase  response.  This  concept  is  covered  in  detail  in 
Chapter  V,  System  Testing  and  Performance. 

The  second  function  accomplished  on  the  synchronizer  board  is  the  detection 
of  the  end  of  the  svncbaud.  After  reference  register  loading,  the  correlator  chip 
is  ready  to  receive  the  syncbaud.  The  TMC2221  provides  a  correlation  sum  from 
zero  to  128  on  an  eight-bit  wide  bus.  By  AN  Ding  together  the  sum  bits,  any 
correlation  score  can  be  detected.  This  is  especially  important  because  rarely  does 
the  system  achieve  a  perfect  128  point  correlation.  Typically,  correlation  sums  fall 
in  the  range  112  to  120.  This  is  due  to  system  noise  and  channel  filtering  per¬ 
turbing  the  waveform  so  that  the  polarities  of  the  received  syncbaud  do  not  cor- 
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respond  exactly  to  the  polarities  of  the  transmitted  syncbaud.  Another  method 
of  obtaining  the  sampling  trigger  is  to  convert  the  correlator  sum  on  the  eight-bit 
sum  bus  through  a  D/A  converter.  The  analog  correlation  values  could  then  be 
compared  to  a  threshold  to  determine  when  to  start  sampling  the  incoming 
waveform.  Only  the  former  method  was  implemented  in  this  thesis  project. 

C.  REAL  TIME  SIGNAL  PROCESSING 

Another  goal  of  this  thesis  was  to  process  the  incoming  bauds  in  real-time. 
For  this  project,  real-time  processing  is  defined  as  being  able  to  complete  the  Fast 
Fourier  Transform  (FFT)  and  the  decoding  of  each  baud  in  one-half  the  time  it 
takes  to  transmit  it.  The  limiting  factor  here  is  the  FFT,  since  the  phase  decoding 
is  a  simple  table  lookup  function  and  occurs  almost  instantaneously  in  compar¬ 
ison  with  the  FFT.  Thus,  only  the  FFT  speed  will  be  discussed  here. 

It  takes  67  milliseconds  to  transmit  a  baud  size  of  4096  samples  when  sam¬ 
pling  at  fx  =  61440  Hz  .  Using  the  Borland  Turbo  Pascal  Toolbox  on  an  IBM 
PC  AT  compatible  machine  running  at  eight  MHz,  with  an  INTEL  80287  nu¬ 
meric  coprocessor,  the  IFFT/'FFT  operations  each  take  one  minute.  Since  up  to 
eight  size  4096  bauds  may  be  received  at  any  time  in  a  packet,  it  could  take  up 
to  eight  minutes  to  recover  the  message  transmitted.  Clearly,  this  is  too  much  time 
to  pay  for  32  kbytes  of  data.  The  key  to  making  MFM  competitive  with  high 
speed  modems  is  to  make  this  operation  occur  in  real-time  with  vector  signal 
processors  (VSPs). 

Several  companies  have  commercially  available  VSPs.  These  chips  are  avail¬ 
able  as  is  or  as  the  featured  component  in  sophisticated  combination  data 
acquisition/signal  processing  boards.  After  investigation  of  several  of  the 
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available  products,  the  PL  1250  plug-in  VSP  board  was  purchased  from 
Eighteen- Eight  Laboratories.  This  VSP  runs  at  12.5  MFLOPS,  using  the  AT&T 
DSP32  chip  as  its  floating  point  engine.  The  board  has  60  kbytes  of  onboard 
static  RAM  for  storing  data  sent  to  it.  The  data  is  transferred  to  the  board  via 
either  register-to-register  transfers  or  DMA.  It  computes  a  size  4096  real  FFT  in 
29  milliseconds,  well  within  the  definition  of  real-time  used  here.  The  board  is 
controlled  using  calls  from  programs  written  in  any  high  level  language.  This  card 
greatly  reduced  the  FFT/'IFFT  processing  time  but  introduced  another  factor 
that  adversely  affected  the  real-time  processing  goal. 

The  data  acquired  and  brought  into  the  receiver  PC  memory  had  to  be  sent 
to  the  VSP.  The  PL1250  calls  initiated  this  process  via  either  register-to-register 
or  DMA  transfers.  In  register-to-register  mode,  the  maximum  data  transfer  oc¬ 
curs  at  500,000  bytes  per  second,  while  DMA  transfers  occur  at  283,000  bytes 
per  second  on  an  AT  class  machine  at  10  MHz  [Ref.  6:  p.5-2].  Thus  in  the 
fastest  of  these  modes,  the  roundtrip  time  for  a  4096  size  baud  (using  32  bit  real 
values  to  represent  each  sample)  is  65  msecs.  This  is  slower  than  the  FFT  com¬ 
putation.  Future  MFM  research  must  address  implementing  the  entire 
acquisition/signal  processing  algorithm  on  a  single  PC  plug-in  card  in  order  to 
compete  with  other  high  speed  modem  techniques. 
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IV.  IMPLEMENTATION 


The  architecture  of  the  MFM  packet  communications  system  requires  an  ex¬ 
tensive  software/hardware  interface.  Some  of  the  required  operations  were  writ¬ 
ten  in  Turbo  Pascal,  while  others  were  coded  in  assembly  language.  It  is 
worthwhile  to  note  that  the  transmitter  uses  an  NPS-designed  and  built  D/A 
converter  board  that  is  driven  with  assembly  language  routines  while  the  receiver 
is  implemented  with  the  Metrabyte,  Inc.,  DASHI6F  data  acquisition  board  and 
Pascal  language  calls  to  accomplish  A/D. 

An  important  difference  between  the  current  transmit  and  receive  devices  is 
that  the  transmitter  is  configured  to  transmit  up  to  61440  byte  samples  while  the 
DASH16F  can  only  acquire  32767  16-bit  samples.  The  actual  Metrabyte  A/D 
conversion  uses  twelve  bits,  while  the  remaining  four  bits  of  the  16-bit  word  are 
used  to  record  the  A/D  channel  used  (up  to  16  channels  available).  This  factor 
comes  into  play  in  the  software  because  the  receiver  can  only  acquire  one-half  of 
the  maximum  length  transmittable  message.  In  this  project,  the  message  was  not 
altered  to  accommodate  the  receiver;  any  data  beyond  the  32767fA  sample  was 
ignored. 

The  programs  discussed  here  are  similar  in  name  and  form  to  those  developed 
by  Gantenbein  [Ref.  1].  They  have  been  modified  for  the  implementation  of  the 
VSP  and  to  encode  and  decode  DI6QAM. 
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A.  TRANSMITTER 


The  transmit  program  is  named  QAMXMIT.  It  uses  either  the  magnitude 
or  Gray-encoded  16QAM  constellations  in  its  encoding  loop.  The  functions  of  the 
program  are  split  into  Pascal  procedures  (subroutines),  which  are  shown  in  Fig¬ 
ure  10.  SYNCBAUD  generates  the  syncbaud.  It  uses  the  same  parameters  for 
amplitude  and  phase  each  time  it  is  called,  so  the  sample  value  polarities  remain 
constant  regardless  of  the  transmitted  message.  Note,  however,  that  samples  of 
the  syncbaud  are  arbitrary.  These  samples  can  be  configured  to  establish  con¬ 
trolled  access  for  the  message  transmission  and  reception. 

SYNCBAUD,  in  turn,  uses  CNVTTOTIME,  which  calls  the  PL1250 
processor  to  perform  the  IFFT  on  the  syncbaud  tone  phase/amplitude  array. 
Next,  SELECTBAUD  assigns  values  to  ku  k2,  and  k„  based  on  the  desired  baud 
length,  such  that  the  signal  will  fall  in  the  16-20  kHz  passband. 
TAILORPACKET,  using  k„  k2 ,  and  k„  from  SELECTBAUD,  determines  the 
number  of  bauds  that  can  be  placed  into  one  segment  of  PC  memory,  64  kbytes. 
It  also  determines  how  much  of  the  message,  stored  in  the  ASCII  file 
"MESSAGE.DAT",  may  be  transmitted  in  the  packet. 

D1FFENCODE  performs  the  encoding  process  as  described  previously,  in¬ 
cluding  the  placement  of  the  complex  conjugates  in  the  second  half  of  the  IFFT 
array.  It  also  calls  CONVTTOT1ME  to  perform  the  IFFT.  Next,  SCALEDATA 
takes  the  time  domain  samples  and  places  them  in  a  vector  of  61440  byte  values 
(the  'broadcast  vector")  to  be  clocked  out  through  the  D/A  converter. 
SCALEDATA  also  checks  each  sample  value  to  ensure  that  it  is  within  the  range 
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Figure  10.  Block  diagram  of  procedures  in  QAMXMIT 


of  the  eight-bit  D/A  converter.  If  the  value  is  out  of  range,  an  error  message  is 
provided  and  the  program  continues  to  run.  The  encode/'IFFT/scaledata  loop  is 
repeated  until  either  MESSAGE.DAT  is  exhausted  or  the  broadcast  vector  is 
filled.  When  the  broadcast  vector  is  ready,  DMAINIT  clocks  the  broadcast  vec¬ 
tor  through  the  D/A  converter.  DMASTOP  is  called  to  halt  the  DMA  process 
after  the  128"1  reference  value  has  been  placed  on  the  PC  data  bus  for  the 
synchronizer  board. 

B.  RECEIVER 

The  receiver  program  is  called  QAMREC.  This  program  acquires,  processes, 
and  decodes  a  received,  magnitude  or  Gray-encoded,  MF16QAM  signal.  Its 
block  diagram  is  shown  in  Figure  11.  Exactly  as  in  the  transmitter,  QAMREC 
begins  by  initializing  the  PL1250  processor.  Next,  GETDMABUFFER  gets  a 
pointer  to  an  array  which  will  be  used  to  store  the  received  samples. 
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Figure  1 1.  Functional  block  diagram  of  program  QAMREC 

PACKETSETUP  assigns  the  kl  and  k2  values  and  determines  the  maximum 
number  of  bauds  recoverable  based  on  kx,  the  receive  baud  length  (which,  of 
course,  must  be  the  same  as  the  transmit  baud  length).  ACQUIREDATA  then 
initializes  the  DASH16F  A/D  board  by  sending  it  set  u*.  parameters  via  h:gh- 
level  language  calls  from  a  software  package  designed  to  drive  the  board  (the 
software  was  developed  by  Quinn-Curtiss  Software,  Inc.).  The  parameter  values 
passed  to  the  board  are  the  board  address,  the  base  storage  address  of  the  re¬ 
ceived  data,  DMA  mode,  and  sampling  trigger  source.  The  DASH16F  stores  the 
received  data  in  unipolar  format,  with  16  bits  per  sample.  The  four  most  signif¬ 
icant  bits  of  each  sample  word  are  used  to  record  the  A/D  channel  used. 

Routine  CONVERTDATA  removes  the  DC  value  introduced  in  the 
DASH16F  sampling  process  and  strips  off  the  four  most  significant  bits.  The  re¬ 
sulting  12  bit  words  are  returned  in  a  real-valued  array  and  the  PL  1250  is  called 
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to  perform  the  FFT.  The  PL  1250  returns  values  in  the  first  half  of  a  complex¬ 
valued  array.  DIFFDECODE  retrieves  the  phase  and  magnitude  difference 
between  adjacent  harmonics  kt  and  ki+]  and  determines  the  ASCII  representation 
for  every  two  symbols  decoded.  The  ASCII  characters  are  displayed  on  the  re¬ 
ceiver  monitor.  The  CONVERTDATA,  FFT,  and  DIFFDECODE  functions 
are  performed  baud  by  baud  until  the  entire  received  message  is  displayed. 

C.  SYNCHRONIZER 

The  synchronizer  board  is  controlled  by  the  software  module,  SYNCLOAD. 
SYNCLOAD  calls  SYNCINIT,  a  slightly  modified  version  of  DMAINIT. 
SYNCLOAD  sets  up  the  parameters  that  will  be  used  in  sending  the  correlation 
reference  values  to  the  TMC2221.  These  parameters  are  the  number  of  bytes  to 
transfer  and  the  location  of  the  bytes  in  memory.  The  real  workhorse  of 
SYNCLOAD  is  the  assembly  language  routine  SYNCINIT,  which  initializes  the 
INTEL  8237  DMA  controller  on  the  PC  motherboard  to  load  the  data  from  PC 
memory  to  the  TMC2221  correlator  chip.  The  PC  bus  lines  driven  by  the  DMA 
controller  are  used  to  control  the  loading  of  the  reference  register  in  the 
TMC2221.  Listings  of  these  programs  are  given  in  Appendices  B  and  C. 
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V.  SYSTEM  TESTING  AND  PERFORMANCE 


The  system  phase  response  and  its  impact  on  differential  phase  encoding  is 
described  in  this  section.  The  performance  measure  signal-to-noise  ratio  out 
(SNReu)  is  defined  and  theoretical  and  observed  values  are  compared.  Last,  the 
number  of  bit  errors  for  10,000  transmitted  bits  for  various  baud  sizes  and  noise 
levels  is  presented. 

A.  PHASE  RESPONSE 

The  phase  responses  of  the  system  for  two  different  cases  of  synchronization 
delay  are  presented  in  Figure  12.  The  synchronization  delay  is  the  number  of 
clock  cycles  between  the  end  of  the  syncbaud  at  the  receiver  and  the  start  of  the 
sampling.  Phase  response  is  defined  as  the  difference  between  the  received  phase 
before  decoding  and  the  transmitted  phase  after  encoding.  It  is  measured  for  each 
tone  in  the  baud,  and  represents  the  amount  of  phase  shift  introduced  in  each 
tone  by  the  channel.  It  is  seen  here  that  for  the  16-20  kHz  passband  (simulated 
by  an  LC  filter),  the  magnitude  of  the  phase  difference  increases  with  frequency 
for  the  zero  delay  case. 

If  different  amounts  of  phase  change  are  introduced  for  each  tone  in  a  baud, 
there  is  a  greater  likelihood  that  the  phase  difference  between  adjacent  tones, 
which  carries  the  information,  will  be  decoded  incorrectly.  If  the  amount  of  phase 
change  added  to  each  tone  by  the  channel  is  identical,  then  the  received  adjacent 
tone  symbol  phase  difference  will  only  be  due  to  the  encoding  process.  The  time 
delay/phase  shift  property  of  Fourier  transforms, 
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PHASE  RESPONSE 


o 


FFT  BIN  NUMBER  IAF-60HZ) 

Figure  12.  Phase  response  for  different  delays 

x(t  -  Td)  ~  X{f)e-jw\  (5) 

shows  that  linear  phase  changes  can  be  controlled  by  adjusting  the  synchroniza¬ 
tion  delay.  For  a  delay  of  one  clock  cycle  in  the  sampling  of  the  incoming 
waveform,  the  phase  shift  at  harmonic  number  290  is  1.8  radians.  This  compares 
well  with  the  phase  difference  between  the  two  curves  in  Figure  12,  1.9  radians 
at  harmonic  number  290.  The  discontinuity  in  the  plot  of  the  zero  delay  case  is 
due  to  "wrap-around"  of  the  phase  difference  from  -n  radians  to  n  radians. 
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The  phase  response  can  be  manipulated  by  selecting  which  128  samples  of  the 
syncbaud  are  used  as  the  correlation  sequence  values  in  the  TMC2221.  Note  that 
the  TMC2221  introduces  a  six  clock-cycle  pipeline  delay  between  the  time 
correlation  occurs  and  the  time  that  the  correlation  sum  is  available  on  the 
eight-bit  correlation  sum  bus  of  the  TMC2221  [Ref.  7J.  Thus  to  obtain  a  syn¬ 
chronization  delay  of  zero  it  is  necessary  to  use  syncbaud  values  123  through  250 
as  the  reference  sequence.  The  maximum  correlation  sum  is  then  output  from  the 
TMC2221  as  the  256rt  sample  of  the  syncbaud  enters  the  TMC2221.  This  sum 
is  then  used  to  trigger  the  A/D  converter  to  begin  sampling  on  the  first  sample 
of  the  first  information  bearing  baud.  From  measurements,  however,  the  best 
phase  response  occurred  when  syncbaud  values  1 24  through  25 1  were  used  as  the 
reference  sequence.  This  produces  an  effective  synchronization  delay  of  one. 
Thus,  the  phase  response  can  be  manipulated  to  optimize  the  channel  for  differ¬ 
ential  encoding  schemes.  A  synchronization  delay  of  one  was  chosen  for  the  test¬ 
ing  described  in  the  remainder  of  this  chapter. 

B.  SIGNAL-TO-NOISE  RATIO  AND  BIT  ERRORS 

The  differentially  encoded  MF16QAM  system  was  tested  for  its  performance 
in  an  AWGN  environment.  Different  levels  of  input  noise  were  introduced  into 
the  system  and  the  SAT?,,,,,  was  calculated  for  each.  The  SAT? euJ  is  defined  as  the 
square  of  the  normalized  mean  level  of  the  complex  multiplied  adjacent  tones  di¬ 
vided  by  their  variance.  The  decoding  process  involves  the  complex  multipli¬ 
cation  of  adjacent  tones.  This  results  in  the  creation  of  one  of  three  tone  product 
levels.  The  three  possible  product  magnitudes  result  when  adjacent  tone 
amplitudes  are:  both  small  (inner  magnitude),  one  small  and  one  large  (middle 


28 


magnitude),  and  both  large  (outer  magnitude).  The  SNR nt  must  be  calculated 
separately  for  each  product  magnitude  level.  Further,  though  the  received  baud 
has  the  same  relative  tone  levels  (\X  and  2X)  as  the  transmitted  baud,  the  actual 
levels  are  imbedded  in  noise.  Since  the  actual  magnitudes  of  these  complex 
multiplications  have  been  corrupted  by  noise,  the  means  and  variances  have  to 
be  computed  conditioned  on  the  actual  magnitudes  of  the  corresponding  adjacent 
transmitted  tones. 

The  bit  error  rate  and  SNR,ut  calculations  are  performed  in  program  SNR, 

which  is  listed  in  Appendix  F.  The  general  steps  of  the  program  are  as  follows: 

•  Export  the  magnitude  and  phase  of  the  recorded  encoded  transmitted  tones 
from  the  transmitter  PC  memory  to  a  File  in  the  receiver  PC  memory. 

•  For  each  baud,  decode  the  received  tones  as  previously  described. 

•  While  the  received  tones  are  being  decoded,  decode  the  transmitted  tones. 

•  Compare  the  decoded  transmitted  and  received  symbols;  if  the  received 
symbol  differs  from  the  transmitted  symbol,  record  the  number  and  type 
(magnitude  or  phase)  of  the  bit  errors. 

•  Combine  the  decoded  symbols  into  ASCII  characters  and  display  them  on 
the  screen  in  different  colors,  depending  on  the  type  of  error  incurred  in  the 
symbol  (if  no  error  occurs,  the  character  is  white). 

•  While  this  joint  decoding  runs,  sort  the  received  complex  multiplied  tones 
into  three  bins,  each  representing  the  one  of  the  three  possible  magnitudes 
resulting  from  multiplying  adjacent  tones.  The  sorting  is  done  based  on  the 
transmitted  values. 

•  After  this  sorting  is  completed,  Find  the  means  and  variances  of  each  of  the 
three  bins  and  calculate  the  SN Rour 

The  curves  in  Figure  13  show  the  SNR out  for  the  inner,  middle,  and  outer 
product  magnitudes  for  a  baud  size  of  1024  with  various  levels  of  SNRn  .  Note 
that  the  SNR,ut  is  lower  as  the  magnitude  product  decreases.  This  is  due  to  the 
AWGN  affecting  each  of  the  product  magnitude  levels  with  equal  measure,  thus 
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causing  the  largest  degradation  when  two  adjacent  transmitted  tones  are  both  in 
level  IX. 

The  SNR„U,  curves  for  the  inner  product  magnitude  level  are  presented  in 
Figure  14.  Also  shown  on  this  curve  is  the  theoretical  SNR for  inner  magni¬ 
tudes  as  given  by  {Ref.  8] 


SNR0Ul  =  SNRJ 


c+  1 


). 


(6) 
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Figure  14.  Output  SNR  for  inner  powers,  all  baud  sizes 


where  c  is  defined  as  the  ratio  of  the  outer  tone  power  to  inner  tone  power  and 
the  SNR  values  are  in  linear  form.  In  Figure  14  it  is  seen  that  the  theoretical 
performance  from  (6)  is  greater  than  the  observed  values  at  high  SNR,„.  This  is 
due  to  the  system  noise  limiting  the  overall  performance  of  the  system.  With  no 
AWGN  introduced,  the  maximum  SNR^,  observed  wras  about  20  dB;  thus 
AWGN  can  only  degrade  the  SNR0U,  below  this  value. 

As  expected,  the  longer  bauds  have  generally  better  performance  (for  5ArR,„ 
above  18  dB).  This  is  because  of  the  smaller  frequency  spacing  between  adjacent 
tones.  At  low  SAR,„,  the  performance  curves  exceeded  the  theoretical.  This  could 
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possibly  be  due  to  positive  correlation  of  the  noise  of  the  adjacent  tones,  which 
would  improve  system  performance  above  the  theory  which  is  based  on  statis¬ 
tically  independent  noise  in  adjacent  frequency  bins  of  the  FFT. 


A  more  quantitative  assessment  of  the  system  performance  for  Gray-encoded 
MFM  is  given  in  Table  2  below,  which  lists  the  number  of  phase  and  magnitude 
bit  errors  for  10,000  transmitted  bits  for  each  baud  size  at  different  SNR»  levels. 
The  larger  the  baud  size  and  the  greater  the  SN the  better  the  system  per¬ 
formance,  as  expected.  Note  that  the  types  of  bit  errors  are  relatively  evenly  split 
between  magnitude  decoding  errors  and  phase  decoding  errors.  Further,  because 
the  phase  errors  are  almost  all  of  the  one  bit  type,  it  is  implied  that  most  of  the 
phase  decoding  errors  were  due  to  phases  being  decoded  one  sector  away  from 
their  actual  transmitted  sectors.  For  this  reason,  the  alternate  (binary)  bit  mag¬ 
nitude  encoding  scheme  was  not  tested,  as  it  would  clearly  introduce  more  bit 
errors. 
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Table  2.  NUMBER  OF  MAGNITUDE  AND  PHASE  DECODING  ERRORS  IN 
10000  TRANSMITTED  BITS  USING  GRAY-ENCODING 


Baud  Size 

256 

512 

1024 

2048 

4096 

SNRIN(dB) 

25 

Mag  bit  errors 

98 

40 

30 

12 

6 

1  bit  phase  errors 

135 

50 

36 

5 

3 

2  bit  phase  errors 

8 

0 

2 

0 

0 

3  bit  phase  errors 

2 

0 

0 

0 

0 

SNRIN(dB) 

20 

Mag  bit  errors 

146 

109 

58 

26 

18 

1  bit  phase  errors 

215 

143 

60 

19 

26 

2  bit  phase  errors 

3 

4 

0 

0 

0 

3  bit  phase  errors 

0 

1 

0 

0 

0 

SNRiN(dB) 

15 

Mag  bit  errors 

261 

222 

222 

201 

193 

1  bit  phase  errors 

376 

275 

283 

191 

191 

2  bit  phase  errors 

5 

2 

2 

0 

0 

3  bit  phase  errors 

0 

0 

0 

0 

0 

SNRIN(dB) 

10 

Mag  bit  errors 

565 

552 

503 

545 

542 

1  bit  phase  errors 

662 

683 

682 

630 

371 

2  bit  phase  errors 

32 

31 

16 

21 

16 

3  bit  phase  errors 

1 

2 

0 

1 

2 

33 


VI.  CONCLUSIONS  AND  FUTHER  STUDY 


This  report  presents  three  advances  which  were  made  in  a  working  MFM 
packet  communication  system.  The  signal  processing  operations  were  made 
real-time  through  the  use  of  vector  signal  processors.  The  synchronization  process 
was  improved  with  the  design  and  implementation  of  a  PC  plug-in  card  that  al¬ 
lows  experimentation  with  the  synchronization  delay  and  software  loading  of  the 
reference  waveform.  A  differential  16  quadrature-amplitude-modulation  encod¬ 
ing  technique  was  developed  and  used  that  doubled  the  data  rate  of  the  MFM 
packet  communications  system  at  NPS.  The  system  was  tested  for  bit  error  rate 
and  signal-to-noise  ratio  performance  in  various  levels  of  AWGN.  Experimental 
results  show  that  MFM  in  its  current  development  can  compete  with  other  high 
speed  modem  techniques  at  input  signal  to  noise  ratios  of  15  dB  and  higher. 

MFM  has  many  areas  in  which  performance  can  be  further  enhanced. 
Larger  signal  constellations  can  be  developed  that  will  increase  the  bandwidth 
efficiency  of  the  MFM  system.  In  concert  with  this,  larger  FFTs  can  be  imple¬ 
mented  to  increase  the  channel  data  volume.  The  D/A  conversion  in  the  trans¬ 
mitter  and  the  data  acquisition  process  in  the  receiver  can  be  expanded  to  12-  or 
16-bit  operations  to  increase  sample  resolution  which  will  undoubtedly  be  neces¬ 
sary  for  the  larger  signal  constellations.  In  the  receiver,  to  achieve  real-time  de¬ 
coding,  the  FFT  operation  should  be  done  baud  by  baud  as  the  data  is  being 
acquired.  While  a  baud  is  being  received,  the  FFT  can  be  performed  on  the 
previously  acquired  baud.  This  is  the  ideal  real-time  decoding  method  and  is 


34 


required  for  a  continuous  transmission  mode.  Finally,  the  receive  and  transmit 
operations  should  be  condensed  onto  one  PC  card  to  enable  each  PC  to  be  a 
transmit  and  a  receive  terminal.  This  would  make  it  possible  to  implement 
networks  with  full  duplex  MFM  modems. 
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APPENDIX  A.  SYNCHRONIZER  DESIGN 


The  synchronizer  has  two  operating  modes,  reference  register  loading  and 
syncbaud  detection.  The  registers  in  the  TMC2221  used  in  this  project  are  as 
follows  [Ref.  7]: 

•  A,  -  128  point  serial  shift  register  into  which  the  data  from  the  hard  limiter 
is  placed.  The  data  is  clocked  in  using  a  61440  Hz  clock.  The  data  in  this 
register  is  compared  with  the  reference  sequence  in  R,  on  each  clock  cycle  to 
determine  the  correlation  score. 

•  B,  -  128  point  serial  shift  register  which  is  used  to  preload  the  reference  syn¬ 
chronization  values.  The  values  are  clocked  in  by  the  61440  Hz  clock  used 
above. 

•  R,  -  reference  latch.  128  parallel  load  register  into  which  the  reference  values 
are  placed.  The  reference  values  are  loaded  in  from  register  B(. 


The  TMC2221  pins  used  are: 

•  LDR  -  load  reference  control  which  copies  the  contents  of  serial  input  regis¬ 
ter  Bf  into  the  reference  register  R,  for  correlation.  If  LDR  was  high  on  the 
previous  clock  cycle,  the  contents  of  input  register  B,-  are  copied  into  input 
register  R(.  to  be  used  as  the  reference  sequence. 

•  A I  -  Data  input  from  the  hard  limiter  is  clocked  in  here  for  comparison  with 
the  reference  sequence. 

•  BY  -  port  where  the  reference  sequence  is  loaded  into  register  B,. 

•  RE1  and  RE2  -  are  used  to  select  either  BY  or  BX  as  the  reference  value 
input  pin.  Both  are  tied  high  to  select  BY  as  the  reference  input  pin. 

•  DM7  through  DM0  -  the  correlation  sum  is  provided  on  this  eight-bit  bus 
after  a  six  clock-cycle  pipeline  delay 

A.  REFERENCE  REGISTER  LOADING 

The  128-point  sequence  of  ones  and  zeros  that  make  up  the  reference  se¬ 
quence  are  stored  in  the  ASCII  file  SYNCVALS.DAT  in  the  receiver  PC.  The 
values  are  read  out  to  the  PC  bus  via  DMA  reads  in  the  single  byte  transfer 
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mode.  In  each  DMA  cycle  the  IOW  (input/output  write)  line  is  drawn  low  when 
the  ASCII  representation  of  each  reference  value  is  placed  on  the  PC  bus  data 
lines  D7-D0.  The  IOW  line  is  used  to  drive  a  latch  (LS374)  to  hold  the  value  of 
the  DO  bit  position  of  the  ASCII  value  which  is  a  one  for  an  ASCII  '1*,  and  is  a 
zero  for  an  ASCII  "0".  This  latched  value  is  then  loaded  into  the  B,  register  of  the 
TMC2221  on  the  next  positive  going  edge  of  the  61440  Hz  clock. 

The  entire  reference  sequence  must  be  placed  into  the  reference  register  R, 
when  the  128,A  value  is  read  out  of  the  file  SYNCVALS.DAT.  This  is  accom¬ 
plished  using  the  DACK3  signal  which  is  driven  low  just  before  IOW  is  driven 
low  on  each  DMA  read.  The  DACK3  signal  drives  the  clock  input  to  a  J-K  flip- 
flop  which  has  been  configured  as  a  toggle  flip-flop  (J  =  K  =  1).  The  output  of  the 
flip-flop  then  changes  state  on  each  DMA  read  and  is  used  to  drive  LDR.  LDR 
is  high  on  every  odd-numbered  DMA  read  and  low  on  every  even-  numbered 
DMA  read,  so  the  128  values  in  B,  are  parallel  loaded  into  R,  on  every  even- 
numbered  DMA  read.  After  the  128,A  DMA  byte  transfer,  DMASTOP  is  called 
to  halt  the  DMA  read  request.  At  this  point,  the  correct  reference  sequence  has 
been  placed  in  R, .  The  schematic  is  given  in  Figure  15. 

B.  CORRELATION 

After  reference  register  loading,  the  synchronizer  board  is  ready  to  perform 
correlation  on  the  hardlimited  signal  from  the  link.  The  correlation  sum  is  pro¬ 
vided  in  the  range  of  zero  to  255  on  DM7  through  DM0.  Synchronization  only 
occurs  at  the  beginning  of  packet  reception,  after  which  any  size  packet  may  be 
continuously  sampled  and  processed  (given  an  MFM  system  that  operates  in 
real-time  as  discussed  in  Chapter  VI).  The  current  synchronizer  board  must  be 
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Figure  15.  Synchronizer  schematic 


reset  after  the  reception  of  each  packet  by  pushing  a  reset  switch  as  shown  in  the 
schematic.  Implementation  of  a  design  for  continuous  reception  of  packets  is  a 
future  research  topic. 
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APPENDIX  B.  SYNCLOAD 

program  sync load; 

(  *****************************  irk** irk*  ***  **  A  *****i *****  *  **  kk  * 

Purpose:  Loads  the  correlation  reference  values  into  the 
TMC2221  correlator  chip  using  DMA  reads  in  the 
single  byte  transfer  mode. 

Inputs:  The  128  correlation  reference  values  must  be  in  the 
file  SYNCVALS.DAT. 

Output:  The  syncvals  (byte  values)  are  placed  on  the  PC  bus 

irk irtrk k** **  *  A*  A-A  * ** ****** *********  *  AA-A  Ai A  jHHHi *************** A1 ) 


uses  crt; 


type 

reference_array  =  array[ 1. . 128]  of  byte; 


var 

j 

reference_values 

num_ref_vals 

vals 

data 

#$L  SYNCINIT 


:  integer; 

:  reference_array; 
:  integer; 

: text; 

: byte; 


procedure  SYNCINIT(var  reference_values:  reference_array; 
num_ref_vals:  integer); 

external; 


begin 


end. 


assign(vals , ’ syncvals.  dat' ); 
reset(vals); 
num_ref_vals  :=  127; 
for  j  :  =  1  to  128  do 
begin 

read( vals, data); 
reference_values{  j]  :=  data; 
writeln(reference_values[ j] ); 

end; 

SYNCINIT( ref erence_values , num_ref_vals ) ; 
close(vals); 


(♦declare  128  values  out*) 
(♦read  each  value  from  *) 
(♦SYNCVALS.DAT  *) 


(*  call  DMA  control*) 
(*  program  *) 
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APPENDIX  C.  SYNCINIT 


code  segment  word  public 
assume  cs:  code 
public  syncinit 

; procedure  syncinit(  BLKADDRESS  :  XMITPOINTER; 

;  BYTECOUNT  :  INTEGER); 

;this  procedure  initializes  dma  channel  3  and  sets  the 
; parameters  to  output  the  array  best  by  passing  the  address 
; of  the  array  start  on  the  stack.  BYTECOUNT  is  the  number 
;of  bytes  to  transfer  and  is  pushed  on  the  stack  by  the 
;  calling  program  (from  Ref.  1,  p.  58) 

dma  equ  0 

dmapage  equ  8 Oh 

syncinit  proc  near 
push  bp 

mov  bp,sp  ;use  bp  to  address  stack 

les  di.dword  ptr[ bp+6] ; move  add  of  best  into  es:  di 

mov  al,4bh  ;  dma  chan  3  single  mode,  read,autoinit 

out  dma+ll,al 

out  dma+12,al  ; reset  first/last  ff 

mov  ax,es  ; calc  high  order  4  bits  of  buffer  area 

mov  cl, 4 
rol  ax, cl 

push  ax  ; save  ax  for  dma  start  addr 

and  al,Ofh 

out  dmapage+2,al  ; store  in  ch  3  dma  page  reg 

pop  ax 

and  a 1, Of Oh 

add  ax,di  ; get  page  offset 

out  dma+6,al  ; output  waveform  buffer  start  addr 

mov  al,ah 
out  dma+6,al 

mov  ax ,  [  bp+4]  ; output  dma  byte  count 

out  dma+7,al 
mov  al,ah 
out  dma+7,al 

mov  al,3  ; unmask  ch  3  to  start 

out  dma+10,al 
pop  bp 

ret  6  ;pop  6  bytes  off  stack  for  addr  of  best 

syncinit  endp 
code  ends 
end 


41 


APPENDIX  D.  QAMXMIT 


program  QAMXMIT; 

Transmits  a  syncbaud  and  message  from  file  'MESSAGE.DAT'. 

The  message  is  differentially  encoded  using  16-QAM.  'MESSAGE.DAT1 
is  a  text  file.  It  should  already  exist  before  using  this  program. 
OUTPUT.DAT  is  used  to  collect  data  for  testing*) 


uses  crt,plrte55; 


const 

FIRST_ELEMENT  =  -28929; 


type 

TNvector  =  array[ 0. . 4095]  of  single; 

TNvectorPtr  =  TNvector; 

BCSTARRAY  =  array] FIRST_ELEMENT. .  32767]  of  byte; 


var 

kx , 

kl,k2,I,w, 

NUMBAUDS .MAXNUMBAUDS , 

BAL'DCOUNT ,  BYTECOUNT , 

SYMBOLCOUNT , MAXNUMCHAR , 

MES S AGE S I ZE , dmachn , 

n2p,bk0psz,bklpsz, 

port ,Aadd,proc,ans3  : integer; 

Baad  : word; 

MAGNITUDE, 

CHARACTERS_PER_BAUD , 

PREV_TONE_MAGNITUDE , PREV_PHASE  : s ingle; 

XREALjXIMAG  :  TNvectorPtr; 

INVERSE  : boolean; 

TEMPBYTE .ERROR  : byte; 

BCST  :  BCSTARRAY; 

BYTEFILE  :  file  of  byte; 

TESTFILE  :text; 

ANSWER, 

NEXTCHAR  : char; 

plbuf  : array] 0. . 768]  of  integer; 


(*$L  dmainit*) 
(*$L  dmastop*) 


(* . *) 

procedure  Cnvttotime; 

(♦computes  inverse  FFT,  returns  values  in  XREAL  *) 
type 

pass  =  array]  0. .  8191]  of  single; 
passptr  =  pass; 
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var 

FVALUES  :passptr; 


begin 


new( FVALUES); 

fillchar( FVALUES  ,sizeof( FVALUES  ),0); 
for  i:=  0  to  kx-1  do 


begin 

FVALUES  [  2*i]  :=  XREAL  [  i] ; 

FVALUES  [ 2*i+l]  :  =  XIMAG  [  i] ; 

end; 

plxfto( FVALUES  ,Aadd,2*kx); 
plwtxf; 

vfieee(Aadd,Aadd,2*kx); 
cifft(Aadd,n2p); 
cereal ( Aadd , Badd , kx ) ; 
vt ieee( B add , B add , kx ) ; 
plwtrn; 

plxffm( Badd, XREAL  ,kx); 
plwtxf; 

dispose( FVALUES); 
end; (*Cnvttotime*) 

(* . *) 

(* . *) 

procedure  SyncBaud; 

(*Process  the  synchronization  baud  and  stores  the  256  point 
time  domain  sequence  at  the  beginning  of  the  packet  storage 
area.  *) 


var 

J,  TEMP  :  integer; 

SYNCDATA  :  byte; 

SYNCMAG  :  single; 

syncvals  : text; 
begin 

ass ign( syncvals, 'syncvals.  datr ); 
rewr ite( syncvals ) ; 
kx:  =256; 
n2p:  =8; 

SYNCMAG:  =  MAGNITUDE; 
f illchar(XREAL  ,sizeof(XREAL 
fillchar(XIMAG  ,sizeof(XIMAG 


XREAL 

XREAL 

XREAL 

XREAL 

XREAL 

XREAL 

XREAL 

XREAL 

XREAL 

XREAL 

XREAL 

XREAL 

XREAL 

XREAL 


[68] 

[69] 

[70] 

[71] 

[72] 

[73] 

[74] 

[75] 

[76] 

[77] 

[78] 

[79] 

[80] 
[81] 


-SYNCMAG; 

-SYNCMAG; 

-SYNCMAG; 

-SYNCMAG; 

SYNCMAG 

SYNCMAG 

-SYNCMAG 

SYNCMAG 

-SYNCMAG 

-SYNCMAG 

SYNCMAG 

-SYNCMAG 

SYNCMAG 

SYNCMAG 


XIMAG 

XIMAG 

XIMAG 

XIMAG 

XIMAG 

XIMAG 

XIMAG 

XIMAG 

XIMAG 

XIMAG 

XIMAG 

XIMAG 

XIMAG 

XIMAG 


)  >0); 
)  >0); 
[68]: 

[69] : 

[70] : 

[71] : 

[72] 

[73] 

[74] 

[75] 

[76] 

[77] 

[78] 

[79] 

[80] 
[81] 


=  -SYNCMAG; 
=  -SYNCMAG; 
=  SYNCMAG; 
=  SYNCMAG; 
:=  -SYNCMAG 
: =  SYNCMAG 
: =  SYNCMAG 
: =  SYNCMAG 
:  =  SYNCMAG 
:=  -SYNCMAG 
: =  -SYNCMAG 
:  =  -SYNCMAG 
: =  SYNCMAG 
:=  -SYNCMAG 
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XREAL  [82]  :=  -SYNCMAG  ;  XIMAG  [  82] :  =  -SYNC MAG  ; 
XREAL  [ 83] :  =  SYNCMAG  ;  XIMAG  [  83] :  =  SYNCMAG  ; 


(★complex  cbnjugate  image*) 
for  J  : =  68  to  83  do 
begin 

XREAL  [  256-J] :  =  XREAL  [J]; 

XIMAG  [ 256-J]  :=-XIMAG  [  J] ; 
end; (*for  J*) 

Cnvttotime;  (*compute  the  256  time  domain  values*) 

for  J  :=  0  to  255  do  (*force  values  to  range  0-255*) 
begin  (*for  d/a  conversion*) 

TEMP:  =round( XREAL  [J]  +  126); 
if  TEMP  <  0  then 
TEMP:  =0; 

SYNCDATA:  =TEMP; 

BCST[ J+FIRST_ELEMENT] :  =SYNCDATA; 
end; (*for  J*) 
close(syncvals); 
end; (*SyncBaud*) 

(* . 

procedure  SelectBaud; 

(*SelectBaud  establishes  kx,  kl,  and  k2,  and  n2p*) 
var 

ANSWER  : integer; 

begin 
kx:  =0; 
repeat 

if  kx  <  0  then  writeln('TRY  AGAIN'); 

writeln( ' What  is  the  length  of  the  bauds  (kx)?'); 

writelnCi.  e.  256,  512,  1024,  2048,  4096'); 

readln( ANSWER); 

case  ANSWER  of 

256:  begin 

kl:=68;  k2:=83;  kx:  =256;  n2p:=8; 
end; 

512:  begin 

kl:  =135;  k2:=166;  kx:  =512;  n2p:=9; 
end; 

1024:  begin 

kl:  =269;  k2:=332;  kx:  =1024;  n2p:  =10; 
end; 

2048:  begin 

kl:  =537;  k2:=664;  kx:  =2048;  n2p:=ll; 
end; 

4096:  begin 

kl:  =1073;  k2:=1328;  kx:  =4096;  n2p:=12; 
end; 

end;  (*case  kx*) 

if  kx  =  0  then  kx  :=  -1; 
until  kx  >  0; 


end;  (*SelectBaud*) 

(* . *) 

procedure  TailorPacket; 

(*TailorPacket  sets  the  maximum  number  of  baud  required  to 
encode  the  message*) 

begin 

MESSAGESIZE:  =  f ilesize(BYTEFILE); 

writeln( 'Message  is  '.MESSAGESIZE,1  bytes.'); 

CHARACTERS_PER_BAUD:  =  (k2-kl)/2; 

MAXNUMCHAR:  =  trunc( 61440. 0/kx  *  CHARACTERS_PER_BAUD ) ; 

if  MESSAGESIZE  >  MAXNUMCHAR  then 
begin 

writeln( 'Message  is  to  large.  The  last 
MESSAGESIZE  -  MAXNUMCHAR, 

'  characters  will  not  be  transmitted.'); 

MESSAGESIZE:  =MAXNUMCHAR; 
end; 

MAXNUMBAUDS:  =trunc(MESSAGESIZE  /  CHARACTERS_PER_BAUD); 

if  frac( MESSAGESIZE  /  CHARACTERS_PER_BAUD)  >0.0  then 
MAXNUMBAUDS:  =MAXNUMBAUDS  +  1; 

repeat 

writeln; 

writeln( 'Enter  number  of  ',kx,'  bauds  to  process. 

MAXNUMBAUDS,'  is  the  maximum.'); 
readln(NUMBAUDS) ; 

until  NUMBAUDS  in  [ 1.  .  MAXNUMBAUDS] ; 
end;  (*TailorPacket*) 

(* . *) 

procedure  DiffEncode; 

C*  DiffEncode  differentially  encodes  the  message  on  a  tone-to-tone 
basis.  BYTEFILE  is  read  from  one  byte  at  a  time.  The  byte 
is  isolated  into  two  4-bit  groups.  Then  the  first  three  bits 
in  each  symbol  of  4  bits  are  used  to  determine  the  phase  shift 
between  tones,  and  the  last  bit  of  the  4  bit  symbol  is  to 
determine  the  magnitude  offset  .  The  encoded  tones  are 
converted  to  rectangular  coordinates  and  are  stored  in  the  arrays 
XREAL  and  XIMAG.  Bytes  partially  encoded  are  carried  over  into  the 
next  baud  by  global  variable  TEMPBYTE  *) 


SHORT  VECTOR, LONG  VECTOR , PHASE SHIFT, 

TONE_MAGNITUDE ,TONE_PHASE , 

PREV_TONE_PHASE , PREV_TONE_MAGNITUDE  :  s ingle; 


DELTAPHI ,DELTAMAG  :byte; 

J  : integer; 


begin 

fillchar(XR£AL  ,sizeof(XREAL  ),0); 
fillchar(XIMAG  ,sizeof(XIMAG  ),0); 
L0NG_ VECTOR  :=  MAGNITUDE; 
SH0RT_VECT0R  :  =  LONG  VECT0R*0. 5; 
PR£V_T0NE  MAGNITUDE  :  =  SHORT  VECTOR; 
PREV.PHASE  :=  22.5; 
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XREAL  [kl]  :  =  SHORT.VECTOR  *  cos( 22.  5*pi/180.  0); 
XIMAG  [kl]  :  =  SHORT.VECTOR  *  sin(22. 5*pi/180. 0); 


writeln(TESTFILE .baudcount , '  ',kl,'  ’ ,PREV_TONE_MAGNITUDE, 

’  \PREV_PHASE); 

if  SYMBOLCOUNT  =  0  then 

r ead( by t  e  f i le , TEMPBYTE ) ; 

for  J:  =  (kl  +1)  to  k2  do 
begin 

SYMBOLCOUNT  :=  SYMBOLCOUNT  +  1; 

(*  seperate  magnitude/phase  bits  *) 

if  frac( SYMBOLCOUNT/ 2)  =0.5  then 
begin 

DELTAPHI  :=  (TEMPBYTE  and  $E0)  shr  5; 

DELTAMAG  :=  (TEMPBYTE  and  $10)  shr  4; 

end; 

if  frac( SYMBOLCOUNT/ 2)  =0.0  then 
begin 

DELTAPHI  :=  (TEMPBYTE  and  $0E)  shr  1; 

DELTAMAG  :=  (TEMPBYTE  and  $01); 
if  NOT  EOF(bytefile)  then 
read(bytef ile .TEMPBYTE) 
else 

TEMPBYTE  :=  $02; 

end; 

(*  differentially  encode  the  last  bit  in  the  four  bit  symbol  *) 

if  PREV_T0NE_MAG.  "UDE  =  SH0RT_VECT0R  then 
begin 

case  DELTAMAG  of 

0:  TONE  MAGNITUDE  :=  SHORT_VECTOR; 

1:  T0NE_MAGN I TUDE  :=  L0NG_ VECTOR; 

end; 

end(*  previous  tone  short  case  *) 
else  (*  PREV_T0NE_MAGN ITUDE  =  L0NG.VECT0R  *) 
begin 

case  DELTAMAG  of 

0:  TONE  MAGNITUDE  :=  LONG  VECTOR; 

1:  TONE.MAGN ITUDE  :=  SHORT.VECTOR; 
end; (*  end  previous  tone  long  case  *) 

end; 

(*  Now  use  the  first  three  bits  in  the  symbol  to  determine  the  amount 
of  phase  rotation  to  the  next  encoded  tone  *) 
if  (ans3=l)  then 
begin 

PHASESHIFT  :  =  DELTAPHI  *  45.  0; 

end 

else 

case  DELTAPHI  of 

0:  PHASESHIFT  :=  0; 

1:  PHASESHIFT  :=  45; 


2:  PHASESHIFT  :=  135; 
3:  PHASESHIFT  :=  90; 

4:  PHASESHIFT  :=  -45; 
5:  PHASESHIFT  :=  -90; 
6:  PHASESHIFT  :=  180; 
7:  PHASESHIFT  :=  -135; 
end; (*case  DELTAPHI  of*) 


(*  Now  assign  the  actual  phase  of  the  tone  being  encoded  which  is  a 
function  of  the  previous  phase,  and  the  phaseshift  *) 

TONE.PHASE  :=  PREV.PHASE  +  PHASESHIFT; 
if  TONE_PHASE  >=  360. 0  then 

TONE.PHASE  :  =  T0NE_PHASE  -  360.  0; 

(*  Now  convert  the  magnitude  and  phase  of  the  tone  to  rectangular 
coordinates  *) 

XREAL  [J]  :=  TONE  MAGNITUDE  *  cos(T0NE_PHASE*pi/180); 

XIMAG  [J]  :=  T0NE_MAGN I TUDE  *  sin(TONE_PHASE*pi/180); 

(*  Save  the  newly  encoded  tone's  magnitude  and  phase  for  the  next 
encoding  iteration  *) 

PREV  TONE  MAGNITUDE  :=  T0NE_MAGN I TUDE ; 

PREV_PHASE  :  =  T0NE.PHASE; 

writeln(TESTFILE,baudcount, '  ',J,'  ' ,T0NE_MAGNITUDE, '  ', 

TONE.PHASE); 

end;  (*  end  of  encoding  process  for  one  tone,  encode  next  tone  *) 

(*  Put  the  complex  conjugate  of  the  encoded  tones  in  the  second  half  of 
the  array  before  computing  the  IFFT  for  this  baud  *) 

for  J:=  kl  to  k2  do 
begin 

XREAL  [  kx  -  J] :  =  XREAL  [  J] ; 

XIMAG  [kx  -  J]  :  =-XIMAG  [  J] ; 

end; 

end; (*DiffEncode*) 

(* . *) 

procedure  ScaleData; 

(*ScaleData  converts  each  real  value  in  XREAL  down  to  a  byte 
and  stores  the  byte  in  the  packet  storage  buffer,  BCST. 

INDEX  establishes  the  location  in  the  buffer  of  each  byte 
in  the  packet.*) 
var 

INDEX, J, TEMP  : integer; 

DATA  : byte; 

begin 

for  J  :  —  0  to  kx- 1  do 
begin 
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if  (XREAL  [J]  >  127)  then 
begin 

writeln( ' IFFT  values  of  of  range  of  d/a  converter'); 
halt; 

end; 

TEMP  :  =  round( XREAL  [J]  +  126); 
if  TEMP  <  0  then 
TEMP  :=  0; 

DATA  :=  TEMP; 

(*256  is  added  to  INDEX  to  start  message  bauds 
after  the  sync  baud*) 

INDEX  :=  J+( BAUDCOUNT- 1 )*kx+FIRST_ELEMENT+256; 

BCST[  INDEX]  :  =  DATA; 

(*  if  baudcount  =  1  then 

writeln(testfile, J: 4, '  ', round( XREAL  [J]));  *) 

end; (*for  J*) 
end;  (*ScaleData*) 

(* . *) 

procedure  Dmastop;  external; 

(*Masks  DMA,  stopping  data  transfer.*) 

(* . *) 

procedure  Dmainit(var  BCST:  BCSTARRAY;  BYTECOUNT:  integer);  external; 
(*Assembly  language  procedure  used  to  initialize  and  unmask 
the  DMA  for  data  transfer.  The  source  code  must  be 
converted  to  a  OBJ  file.*) 

(* . *) 

(* . *) 

begin 

dmachn:  =0; 

p 1 in i t ( dmachn , p lbu  f,sizeof(p lbuf ) ) ; 

plslibCC:  PL1250  PLLIB.  15'); 

proc:  =1; 

port: =$0318; 

bkOpsz: =0; 

bklpsz: =1024; 

plspr c( proc, port , bkOpsz ,bklpsz); 

Aadd: =$0400; 

Badd: =$8400; 

new( XREAL); 
new(XIMAG); 

(♦contains  hex  values  to  be  encoded  and  transmitted*) 
assign(BYTEFILE, 'MESSAGE. DAT1 ); 
reset(BYTEFILE); 

(♦Output  file  of  encoded  symbols.  Used  for  system  testing*) 
assign(TESTFILE , ' XMITDAT.  DAT' ); 
rewrite(TESTFILE); 

repeat 

writeln( 'enter  the  type  of  encoding  ( l=magnitude,2=gray) * ); 
readln(ans3); 

writeln( 'Enter  magnitude  of  tones 
(greater  than  65,  less  than  1501)'); 
readln( MAGNITUDE); 
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until  MAGNITUDE  >  0.0; 
writelnC 'Loading  sync  baud. '); 

SyncBaud; 

SelectBaud; 

TailorPacket; 

SYMBOLCOUNT:  =0; 

TEMPBYTE:  =$00; 

writelnC 'Number  of  bauds  is  ', numbauds); 


for  baudcount  :  =  1  to  numbauds  do 
begin 

DiffEncode; 

writelnC 'Performing  IFFT  '.BAUDCOUNT,'  ’ 
NUMBAUDS -BAUDCOUNT,'  left'); 
Cnvttotime; 

ScaleData; 

end;  (*for  BAUDCOUNT*) 


BYTECOUNT  :=  256  +  NUMBAUDS*kx  -  1; 
wr ite ln( bytecount ) ; 

repeat 

writeln( ' Press  return  to  transmit' );  readln; 

DmainitC  BCST , BYTECOUNT) ; 

repeat 

writelnC 'Transmit  some  more?  C*yes  or  no*)  '); 
readlnC ANSWER); 

until  ANSWER  in  [ 'n' , 'N' , 'y' , ' Y' ] ; 

Dmastop; 

until  ANSWER  in  [ ’n' ,'N’] ; 

disposeCXREAL); 

disposeCXIMAG); 

closeCBYTEFILE); 

closeCTESTFILE); 

end. 
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APPENDIX  E.  QAMREC 


program  QAMREC; 

(****rt***rtft*rt»*****^**********ili**iHr  ******  A*  *******  A***********  ****** 

Purpose:  This  program  acquires  the  signal , stores  it  in  a  memory  buffer, 
performs  the  FFTs ,  and  decodes  and  displays  the  received 
symbols  in  the  form  of  ASCII  text. 


Inputs  :  The  inputs  are  taken  form  the  keyboard  in  the  form  of  user 

responses.  The  inputs  are:  (1)  The  number  of  bauds  to  be 

processed,  (2)  the  size  of  the  bauds,  (3)  the  type  of  decoding 

scheme  desired  (magnitude  or  gray). 


Outputs:  The  outputs  are: 

File  RECDAT.  DAT  -  the  real  and  imaginary  parts  of  the  received 
tones 

The  decoded  symbols,  displayed  on  the  receive  CRT. 


uses  Graph,  Crt,  tp55d!6,plrte55; 


(*$I-*) 

(*$R-*) 

const  Max_Buffer  =  65500; 
type 

(*TYPE  for  real  and  imaginary  data  for  FFT  routing*) 
TNvector  =  array[ 0.  . 4095]  of  single; 

TNvectorPtr  =  TNvector;  (*Pointer  for  FFT  data  array 

which  allows  dynamic 
allocation  of  memory*) 


var 

XREAL,  XIMAG 

ERROR,  TEMPBYTE 

J,I ,xradd,xroadd,proc,port, 

kl , k2 , kx .ANSWER ,ERR_COOE , 

BAUDCOUNT , SYMBOLCOUNT , n2p , 

NUM3AUDS , MAXNUMBAUDS , dmachn , 

bkOpsz ,bklpsz , ans2 

MAGNITUDE, PHASE 

DATA VECTOR 

DMAPOINTER 

OUTFILE , recdat 

plbuf 


: TNvectorPtr; 
:  byte; 


: integer; 

:  real; 

:  integer; 

: pointer; 

:  TEXT; 

: array[ 0. . 4095]  of  integer; 


(* . *) 

procedure  PacketSetUp; 

(*  Defines  the  baud  parameters  given  the  desired  baud  length  *) 

begin 

repeat 

clrscr; 


50 


if  kx  <  0  then  writeln('TRY  AGAIN’); 
writeln( 'Enter  baud  size  '); 
readln( ANSWER); 
case  ANSWER  of 
256:  begin 

kx: *  256; 
n2p:  *8; 

end; 

512:  begin 

kx:*  512; 
n2p:  =9; 

end; 

1024:  begin 

kx: *1024; 
n2p:  =10; 

end; 

2048:  begin 

kx: *2048; 
n2p:  =11; 

end; 

4096:  begin 

kx: =4096; 
n2p:  =12; 

end; 

end; (*case*) 

if  kx  =  0  then  kx  :=  -1; 
until  kx  >  0; 

MAXNUMBAUDS  :=  trunc((MAX_BUFFER/2)/kx); 

repeat 

writeln; 

writeln( 'Enter  number  of  ',kx,'  bauds  to  process. 

MAXNUMBAUDS,'  is  the  maximum.  ' ) ; 
readln(NUMBAUDS); 

until  NUMBAUDS  in  [ 1.  . MAXNUMBAUDS] ; 

kl  :=  round(kx  *  67.  0  /  256.0  +  1); 
k2  :=  roundCkx  *  83.  0  /  256.0); 
end;  (*PacketSetUp*) 

(* . *) 

procedure  AcquireData; 

(*AcquireData  initializes  Metrabyte  DASH-16F  data  acquisition 
board,  using  TT00LS  procedure  D16_int  and  D16_ainm.  Data 
transfer  is  controlled  by  the  DMA  controller  and  initialized 
by  Dl6_ainm  and  disabled  by  D16_dma_int_disable.  TTOOLS 
procedures  are  external  procedures  included  by  'uses' 
tp4dl6.  *) 

var 

RATE:  real; 

I,CNT_NUM,  MODE,  CYCLE,  TRIGGER, 

BASE.ADR,  INT.LEVEL,  DMA_ LEVEL, 

B0ARD_NUM ,  CHANLO, 
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OP_TYPE,  STATUS,  NEXT.CNT,  ERR_CODE_S  :  integer; 

begin 

BOARD  NUM  :=  0;  INT.LEVEL  :  =  7;  DMA.LEVEL  :=  1; 

BASE_ADR  :  =  §300; 

D16  in it (BOARD  NUM , BASE_ADR , INT_LEVEL , DMA_LEVEL , ERR_CODE ) ; 
if  £RR_C0DE  <>  0  then 

D16_pr int_error ( ERR^CODE) ; 

CHANLO  : =  0; 

CYCLE:  =0;  (*0-one  sweep  of  the  DMA  1-autoinitialize*) 

TRIGGER:  *0;  (*0  -  external  1  -  internal*) 

CNT_NUM:  =32767;  (*#  of  samples*) 

RATE  :=  10000. 0; (*used  for  internal  trigger*) 

MODE  :=  2;  (*DMA  mode*) 
writeln( 'Ready  to  acquire'); 

Dl6_ainm( B0ARD_NUM , CHANLO , MODE , CYCLE , TRIGGER , CNT_NUM , 

RATE,  DATAVECT0R  ,ERR_C0DE); 
if  ERR_C0DE  <>  0  then 

D 1 6_pr int_er ror ( ERR_C0DE ) ; 

STATUS  :=  11; 

(♦status  indicates  the  progress  of  acquisition.  When  all 
samples  have  been  acquired  status=0*) 
repeat 

D 1 6_dma_ int_s  t at us ( B0ARD_NUM , 0P_TYPE , STATUS ,NEXT_CNT, 

ERR_C0DE_S); 
if  ERR_C0DE_S  <>  0  then 
D 1 6_print_e r  r or ( ERR_C0DE_S ) ; 

until  STATUS  =  0; 
writeln( 'Data  received'); 

if  ERR_C0DE  <>  0  then 

D 16_pr int_error ( ERR_C0DE ) ; 

D16_dma_int_disable(B0ARD_NUM,ERR  CODE); 
if  ERR_C0DE  <>  0  then 

D16_print_error(ERR_C0DE); 

end; (*Acquire*) 

(* . *) 

procedure  ConvertData; 

(*ConvertData  seperates  channel  and  acquired  data,  and  converts  the 
twelve  bit  acquired  data  into  turbo  pascal  six  byte  real  values.  These 
sample  values  are  returned  in  the  array  XREAL  []*) 

var 

AD_DATA:  array[ 0. .  4095]  of  integer; 

I , CHAN_DATA , ERR  CODE, 

SEGMENT? ART , OFFSETPART: integer; 

NEWDATAVECTOR  :  integer; 

TEMPPOINTER  : pointer; 
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begin 

fillchar(xreal  ,sizeof(xreal  ),0); 
fillchar(ximag  ,sizeof(ximag  ),0); 

SEGMENTPART:  =s eg( DATAVECTOR  ) ; 

OFFSETPART:  =ofs(DATAVECTOR  )  +  2  *  kx  *  (BAUDCOUNT  -  1); 
TEMPPOINTER:  =ptr( SEGMENTPART, OFFSETPART); 

NEWDATA VECTOR  :=  TEMPPOINTER; 

dl6_convert  data( 2048 , kx , NEWDATA VECTOR  ,AD  DATA[ 0] , 

CHAN.DATA , 0 , ERR^CODE ) ; 
if  ERR^CODE  <>  0  then 
D16_print_error(ERR_C0DE); 
for  I:  =  0  to  (kx  -  1)  do 
begin 

xreal  [  i]  :  =  AD_DATA[  i] ; 

end; 

end; (*ConvertData*) 

(* . *) 

procedure  DiffDecode; 

(*DiffDecode  differentially  decodes  complex  frequency  domain 
arrays  XREAL  and  XIMAG.  Two  decoded  symbols  are  recombined 
into  a  byte  and  transferred  to  the  screen*) 

var 

I  : integer; 

TEMPREAL .TEMPIMAG , OLDMAG .NEWMAG  :  single; 

B ITS , PHASEB ITS , MAGB IT  :  byte; 

TEMPCHAR  :  char; 

begin 

for  I : =  kl  to  (k2-l)  do 
begin 

SYMBOLCOUNT:  =  SYMBOLCOUNT  +  1; 

(*  save  the  current  and  next  magnitudes  for  future  decoding  *) 

OLDMAG  :=  sqrt(sqr(XREAL  [  I] )  +  sqr( XIMAG  [I])); 

NEWMAG  :=  sqrt(sqr( XREAL  [ 1+1] )  +  sqr(XIMAG  [1+1])); 

(*  complex  multiply  adjacent  tones  to  get  phase  differential  *) 

TEMPREAL  :  =  XREAL  [ I]  *  XREAL  [ 1+1]  + 

XIMAG  [I]  *  XIMAG  [  1+1] ; 

TEMPIMAG  :  =  XREAL  [  I]  *  XIMAG  [ 1+1]  - 
XREAL  [ 1+1]  *  XIMAG  [  I]  ; 

(*  now  rotate  phase  by  22.5  degrees  to  line  up  with  constellation 
phase  sectors  *) 

XREAL  [I]  :=  0.92  *  TEMPREAL  -  0.  38  *  TEMPIMAG; 
XIMAG  [I]  :=  0.92  *  TEMPIMAG  +  0. 38  *  TEMPREAL; 

(*  decode  the  phase  difference  into  the  first  three  bits  of  the  symbol 
to  be  recovered  *) 

if  (ans2  =  1)  then  (*magnitude  encoding  scheme*) 
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begin 

PHASEBITS  :=  $00; 

if  (XREAL  [I]  >0)  and  (XIMAG  [I]  >0)  then 
if  XREAL  [  I]  >  XIMAG  { I]  then 
PHASEBITS  :=  $00 
else  PHASEBITS  :  =  $02; 

if  (XREAL  £  I]  <  0)  and  (XIMAG  [I]  >0)  then 
if  abs( XREAL  [  I] )  >  XIMAG  [I]  then 
PHASEBITS  :=  $06 
else  PHASEBITS  :=  $04; 

if  (XREAL  [I]  <  0)  and  (XIMAG  [I]  <  0)  then 
if  abs( XREAL  [  I] )  >  abs( XIMAG  [I])  then 
PHASEBITS  :=  $08 
else  PHASEBITS  :  =  $0A; 

if  (XREAL  [I]  >0)  and  (XIMAG  [I]  <  0)  then 
if  XREAL  [I]  >  abs( XIMAG  [I])  then 
PHASEBITS  :=  $0E 
else  PHASEBITS  :=  $0C; 

end 

else  (*if  phase  encoding  scheme*) 
begin 

PHASEBITS  :=  $00; 

if  (XREAL  [I]  >0)  and  (XIMAG  [I]  >0)  then 
if  XREAL  1 1]  >  XIMAG  [ I]  then 
PHASEBITS  :=  $00 
else  PHASEBITS  ;=  $02; 

if  (XREAL  [I]  <  0)  and  (XIMAG  [I]  >0)  then 
if  abs ( XREAL  [  I] )  >  XIMAG  [  I]  then 
PHASEBITS  :=  $04 
else  PHASEBITS  :=  $06; 

if  (XREAL  [I]  <  0)  and  (XIMAG  [I]  <  0)  then 
if  abs (XREAL  [ I] )  >  abs (XIMAG  [I])  then 
PHASEBITS  :=  $0C 
else  PHASEBITS  :=  $0E; 

if  (XREAL  [I]  >0)  and  (XIMAG  [I]  <  0)  then 
if  XREAL  [I]  >  abs(XIMAG  [I])  then 
PHASEBITS  :=  $08 
else  PHASEBITS  $0A; 


(*  now  differentially  decode  the  magnitudes  of  the  tones  to  get  the 
fourth  and  last  bit  in  the  symbol  *) 

if  (NEWMAG  >  1.  5*0LDMAG)  or  (NEWMAG  <  2*0IDMAG/3) 
then  MAGBIT  :=  $01 
else  MAGBIT  :=  $00; 


(*  now  jam  all  the  bits  together  *) 
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(*fill  TEMPBYTE  with  two  symbols*) 
if  f rac( SYMBOLCOUNT  /  2)  *  0.  5  then 

TEMPBYTE  :=  ((PHASEBITS  or  MAGBIT)  shl  4); 
if  (frac( SYMBOLCOUNT  /  2)  =  0.0)  then 
begin 

TEMPBYTE  : =  (PHASEBITS  or  MAGBIT)  or  TEMPBYTE; 

TEMPCHAR  :  =  chr( TEMPBYTE); 

write(TEMPCHAR);  (*  put  ascii  character  to  screen*) 

TEMPBYTE:  =0; 
end;  (*if  frac*) 
end;  (*for  I*) 

end;  (*DiffDecode*) 

(* . *) 

begin  (*main  body*) 

(*  initialize  the  PL  processor  *) 
dmachn:  =0; 

plinit(dmachn,plbuf ,sizeof(plbuf )); 

plslib('c:  pl850  pllib. 15'); 

proc:  =  1; 

port: =  $0318; 

bkOpsz: =0; 

bklpsz:  =1024; 

plsprc( proc, port, bkOpsz .bklpsz); 

GetDMABuf  f er ( MAX_BUFFER , DMAPOINTER , ERR_C0DE ) ; 

DATAVECTOR  :=  DMAPOINTER;  (*This  statement  assigns  a 
generic  pointer  to  a  variable  of  a  specific  pointer 
type,  i. e.  integer,  so  that  the  pointer  can  be 
passed  to  the  dl6_ainm  routine.*) 

assign(recdat, 'recdat.  dat' ); 
reu’riteCrecdat); 

writeln( 'Which  decoding  scheme  ( l=magnitude,2=gray) ' ); 
read(ans2); 

new(XREAL); 

nev(XIMAG); 

ERROR  :=  0; 
kx:  =0; 

PacketSetUp; 

SYMBOLCOUNT:  =0; 

TEMPBYTE:  =0; 

AcquireData;  (*AcquireData  samples  input  analog  signal*) 

(*  set  up  address  in  PL  FPP  memory  that  is  the  starting  point  for  where 
the  acquired  data  is  sent  to  for  FFT  processing  *) 
xradd: =$0400; 
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(★begin  baud  by  baud  conversion*) 
for  8AUDC0UNT  :=  1  to  NUMBAUDS  do 
begin 

ConvertData; 

plxfto(xreal  ,xradd,kx); 
plwtxf; 

vfieee(xradd,xradd,kx); 

plwtrn; 

rfft(xradd,n2p); 

plwtrn; 

vt ieee( xr add , xradd , kx) ; 
plwtrn; 

plxffo(xradd, xreal  ,kx); 
plwtxf; 

for  j:=  0  to  kx  div  2  do 
begin 

xreal  [ j]  :=  xreal  [ 2*j] ; 
ximag  [  j]  :  =  xreal  [  2*  j+l] ; 

end; 

ximag  [0]  :  =  0; 
for  i  :  *  kl  to  k2  do 

begin  ,  , 

writeln(recdat,baudcount,  ,1,  , 

XREAL  [ I] : 10: 4, '  ' , XIMAG  [ I] : 10: 4); 

end; 

Diff Decode; 
end; 

close(recdat); 

dispose(XREAL); 

dispose(XIMAG); 

Fr e  eDMABu  f  f  e  r ( MAX_BUFFER , DMAPO INTER , ERR_C0DE ) ; 
if  ERR.CODE  <>  0  then 

D16_print_error(ERR_C0DE); 

end. 


56 


APPENDIX  F.  SNR 


program  snr; 

(  **ft»*A**ftft**»ft***ft******'*»**  A*******  A  ******  A  A'AtM  *******  A*  ******  ********* 

Purpose:  This  program  is  used  to  calculate  the  statistics  of  MFM  for  any 
any  size  baud.  It  determines  the  SNR  out  and  counts  the  number 
and  type  of  bit  errors  in  the  decoding  process. 

Inputs  :  This  program  requires  the  following  data  files  : 

RECDAT.DAT  -  the  received,  decoded  tone  real  and  imaginary  parts 
XMITDAT.DAT  -  the  encoded,  transmitted  tone  magnitudes  &  phases 

Outputs:  The  statistics  and  number  and  type  of  biterrs  are  stored  in  the 
file  output.dat.  In  addition,  the  received  ASCII  characters  are 
displayed  on  the  screen.  If  the  received  ASCII  character  is 
correct,  it  is  displayed  in  white.  Otherwise,  it  is  displayed 
in  one  of  the  following  colors: 

yellow  -  if  a  magnitude  decoding  error  occurred 

green  -  if  a  phase  decoding  error  occurred 

red  -  if  a  phase  AND  magnitude  decoding  error  occurred 

uses  crt,  graph; 


var 

answer ,answer2  :  char; 

i , j , n , rbaud , xbaud , rtone , xtone , dtot , 
baudcount , numbauds , k 1 , k2 , kx , count , 
symbolcount , sector , b , c , d , btot , ctot , 

badbaud.numbits ,badbaud2 ,bj , cj ,dj  : integer; 

xoldmag, roldmag.xtempreal , rtempreal , 

xnewmag , rnewmag , xtempimag , rtempimag , 

totphaserrs , totmagerrs , symerrs , 

smallmag,big,sml ,del ,meanbig,xmagbig,bigmag, 

mbig,msml,obig,osml,mmeanbig,mmeansml , 

omeanbig,oraeansml,mde’.,odel,msnravg,osnravg, 

meansml ,varx,snrbig,sarsml ,xmagsml ,snravg, 

mvarx,ovarx,msnrbig,msnrsml ,osnrbig,osnrsml , 

mxmagbig ,mxmagsml , oxmagbig , oxmagsml , 

totsnr .mtotsnr .ototsnr  : real; 

xbits .xphasebits ,xmagbit ,xtempbyte, 

rbits , rphasebits , rmagbit , rtempbyte , 

phasebitdif f ,magbitdif  f ,hue ,pbdl , 

pbd2 ,pbd3  :  byte; 

xtempchar ,rtempchar  : char; 

xmitdat.recdat, output  : text; 

xreal,ximag,rreal,rimag,xmag,xphase  : array [ 1. . 256]  of  real; 

statmat  :  array]  1. .  8 , 1. .  3]  of 

real; 

recdata  : array]  1. .  48,1. . 100] 

of  real; 

snrin  : string] 4]; 
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(* . *) 

procedure  sort; 

( ■kirk*  A*************** *  ******  ft  *'*  A  **********  A  "A  ******* A A* A* ******* ft Aft****** 

Purpose:  This  procedure  is  used  to  sort  the  complex  multiplied  tones  into 
three  different  magnitude  bins,  inner,  middle,  and  small. 

Inputs  :  The  inputs  are  the  adjacent  transmitted  magnitudes  xmag[ i]  and 
xmagf  i+1] . 

Outputs:  The  output  is  the  global  array  recdata. 

**i^*«*ft ************** ft Aft **************** ************ Aft* Aft* Aft* ft******) 

begin 


if  (xmagf I]  =smallmag)  and  (xmag[  1+1]  =smallmag)  then 
begin 

statmatf  sector ,  1] :  =statmat[  sector,  1]  +1; 
b  :=  round(statmat[ sector, 1] ); 
recdata[ (2*sector)-l,b]  :=RREALf  I] ; 
recdataf (2*sector) ,b] :  =RIMAG[ I] ; 

end 

else  if  (((xmagf I]=smallmag)  and  (xmag[ 1+1]  =bigmag))  or 

((xmagf  I]  =bigmag)  and  (xmagf 1+1] =smallmag)))  then 
begin 

statmat[  sector , 2] :  =statmat[ sector , 2] +1; 
b  :=  round( statmatf sector ,2] ); 
recdata[  (2*sector)-l+16,b] :  =RREAL[  I] ; 
recdataf (2*sector)+16,b]  :=RIMAG[ I] ; 

end 

else  (*  both  xmag[ I]  and  xmagf  I+l]  are  large  *) 
begin 

statmatf  sector ,3] : =statmat[ sector ,3] +1; 
b  :=  round(statmat[ sector, 3] ); 
recdataf (2*sector)-l+32,b]  :=RR£AL[ I] ; 
recdataf ( 2*sector )+32 , bj : =RIMAGf I] ; 

end; 

end; 

(* . 


*) 


begin  (*main  body*) 
clrscr; 

ass ign( output, 'output,  dat' ); 
rewr ite( output ) ; 
assign(xmitdat , 'xmitdat. dat' ); 
reset(xmitdat); 
assign(recdat, 'recdat.  dat' ); 
reset(recdat); 

writeln( 'Enter  the  input  snr'); 
readln(snrin); 

writeln( 'Enter  the  baud  length  '); 
readln(kx); 

writeln( output, 'The  baud  length  is  ',kx,'  and  the  SNRIN  *' ,snrin); 
writeln( 'Enter  the  number  of  bauds  to  be  processed'); 
readln( numbauds ) ; 
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writelnC 'Enter  the  magnitude  of  the  xmit  short  tones'); 

readlnC  smal Imag) ; 

bigmag  :  =  2*smallmag; 

writelnC 'Throw  out  any  bauds  ?  '); 

readln( answer); 

badbaud  : =0; 

badbaud2:  =0; 

if  answer  in  {'y'.'Y'J  then 
begin 

writeln( 'Which  baud  ?'); 
readlnC  badbaud) ; 
writelnC'Any  others  ?'); 
readlnC answer2); 
if  answer2  in  [ ' y ' , ' Y ' ]  then 
begin 

writelnC  'Enter  baud  //'); 

readlnC  badbaud2 ) ; 

end; 

end; 

C*  set  up  for  the  desired  baud  size  *) 
case  kx  of 

256:  begin 

kl:  =68;  k2:  =83; 

end; 

512:  begin 

kl:  =135;  k2:  =166; 

end; 

1024:  begin 

kl:  =269;  k2: =332; 

end; 

2048:  begin 

kl:  =537;  k2:  =664; 

end; 

4096:  begin 

kl:  =1073;  k2:  =1328; 

end; 

end;  C*case  Kx*) 

C*  initialize  overall  statistical  variables  *) 

TOTPHASERRS  :=0; 

TOTMAGERRS  :  =0;  SYMB0LC0UNT  :  =0; numbits:  =0; 
pbdl: =0; pbd2: =0; pbd3  : =0; bj: =0; cj: =0; dj: =0; 

totsnr:  =0;  mtotsnr:  =0; ototsnr: =0; 

(*  Now,  count  bit  errors  baud  by  baud  *) 

C*  read  in  transmit  and  receive  values  *) 
for  j:=  1  to  numbauds  do 
begin 

del:  =0;  big:  =0;  sml:  =0;  meanbig:  =0;  meansml:  =0;  btot:  =0; 
mdel:  =0;  mbig:  =0;  msml:  =0;  mmeanbig:  =0;  mmeansml:  =0;  ctot:  =0 
odel:  =0;  obig:  =0;  osml:  =0;  omeanbig:  =0;  omeansml:  =0;  dtot:  =0 
f illcharCstatmat ,sizeof(statmat) ,0); 
f illcharCrecdata,sizeof(recdata) ,0); 
for  i:=  1  to  k2-kl+l  do 
begin 
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readln(xmitdat, xbaud, xtone,xmag[  i]  ,xphase[  i]  ); 
reading recdat,rbaud,rtone,rreal[ ij ,rimag[ i] ); 

if  (xbaud  <>  rbaud)  or  (xtone  <>  rtone)  then 
begin 

writeln( 'RECDAT  and  XMITDAT  do  not  match'); 
halt; 

end;  (*if  xbaud*) 

(*  convert  the  xmit  vals  to  rectangular  coordinates*) 
xreal[  i]  :=  xmag[  i]  *cos(xphase[  i]  *pi/180); 
ximag[  ij  :  =  xmag[  ij  *sin(xphase[  ij  *pi/180); 
end;  (*for  read  data  files*) 
writeln; 
write( j , '  '); 

(*  Now  commence  conditional  decoding  of  the  received  symbols  *) 
for  I:=  1  to  k2-kl  do 
begin 

SYMBOLCOUNT: =  SYMBOLCOUNT  +  1; 

(*  save  the  current  and  next  magnitudes  for  future  decoding  *) 

XOLDMAG  :=  sqrt( sqr(XREAL[ I] )  +  sqr(XIMAG[  I] )); 
XNEWMAG  :=  sqrt(  sqr(XREAL[  1+1]  )  +  sqr(X!MAG[  1+1]  )  ); 
ROLDMAG  :=  sqrt( sqr(RREAL[  I]  )  +  sqr(RIMAG[ I] )); 
RNEWMAG  :=  sqrt( sqr(RREAL[  1+1]  )  +  sqr(RIMAG[  1+1]  )); 

(*  complex  multiply  adjacent  tones  to  get  phase  differential  *) 

XTEMPREAL  :  =  XREAL[  I]  *  XREAL[ 1+1]  + 

XIMAG[  I]  *  XIMAG[  1+1]  ; 

XTEMPIMAG  :=  XREAL[  I]  *  XIMAG[  1+1]  - 
XREAL[  1+1]  *  XIMAG[  I]  ; 

RTEMPREAL  :=  RREAL[  I]  *  RREAL[  1+1]  + 

RIMAG[  I]  *  RIMAG[ 1+1]  ; 

RTEMPIMAG  :  =  RREAL[ IJ  *  RIMAG[ 1+1]  - 
RREAL[  1+1]  *  RIMAG[  I]  ; 

(*  now  rotate  phase  by  22.5  degrees  to  line  up  with  constellation 
phase  sectors  *) 

XREAL[ I]  :=  0.92  *  XTEMPREAL  -  0. 38  *  XTEMPIMAG; 

XIMAGJ  I]  :=  0.92  *  XTEMPIMAG  +  0.  38  *  XTEMPREAL; 

RREAL[  I]  :=  0.92  *  RTEMPREAL  -  0.  38  *  RTEMPIMAG; 

RIMAGJIJ  :=  0.92  *  RTEMPIMAG  +  0.  38  *  RTEMPREAL; 

(*  decode  the  transmit  phase  difference  into  the  first  three  bits  of  the 
symbol  to  be  recovered  *) 

XPHASEBITS  :=  $00; 

if  (XREALf I]  >  0)  and  (XIMAG[ I]  >  0)  then 
if  XREAL[ I]  >  XIMAG[ I]  then 
begin 

XPHASEBITS  :=  $00; 
sector  :=1; 
sort; 

end 

else 

begin 

XPHASEBITS  :=  $02; 
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sector  :=2; 
sort; 

end; 

if  (XREAL[ I]  <  0)  and  (XIMAG[ I]  >  0)  then 
if  abs(XREAL[ I] )  >  XIMAG[ I]  then 
begin 

XPHASEBITS  :=  $04; 
sector  :  =  4; 
sort; 

end 

else 

begin 

XPHASEBITS  :=  $06; 
sector  :=  3; 
sort; 

end; 

if  (XREAL[ I]  <  0)  and  (XIMAGfl]  <  0)  then 
if  abs(XREAL[ I] )  >  abs(XIMAG[ I] )  then 
begin 

XPHASEBITS  :=  $0C; 
sector  :=  5; 
sort; 

end 

else 

begin 

XPHASEBITS  :=  $0E; 
sector:  =6; 
sort; 

end; 

if  (XREAL[ I]  >  0}  and  (XIMAGfl]  <  0)  then 
if  XREAL[ I]  >  abs(XIMAG[ I] )  then 
begin 

XPHASEBITS  :=  $08; 
sector  :  =  8; 
sort; 

end 

else 

begin 

XPHASEBITS  :=  $0A; 
sector  :=  7; 
sort; 

end; 

(*  decode  the  received  phase  difference  into  the  first  three  bits  of  the 
symbol  to  be  recovered  *) 

RPHASEBITS  :=  $00; 

if  (RREAL[ I]  >  0)  and  (RIMAG[ I]  >  0)  then 
if  RREAL[ I]  >  RIMAG[ I]  then 
RPHASEBITS  :«  $00 
else  RPHASEBITS  :=  $02; 

if  (RREAL[ I]  <  0)  and  (RIMAG[ I]  >  0)  then 
if  abs(RREAL[  I] )  >  RIMAG[ I]  then 
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RPHASEBITS  :=  $04 
else  RPHASEBITS  :=  $06; 

if  (RREAL[ I]  <  0)  and  (RIMAG[ I]  <  0)  then 
if  abs(RREAL[ I] )  >  abs(RIMAG[  I] )  then 
RPHASEBITS  :=  $0C 
else  RPHASEBITS  :=  $0E; 

if  (RREAL[  I)  >  0)  and  (RIMAGf I]  <  0)  then 
if  RREAL[  I]  >  abs(RIMAG[  I]  )  then 
RPHASEBITS  :  =  $08 
else  RPHASEBITS  :=  $0A; 


(*  determine  the  number  of  bit  differences  between  the  received  decoded 
phasebits  and  the  decoded  transmitted  phasebits  *) 

PHASE BITD IFF  :=  XPHASEBITS  xor  RPHASEBITS; 
if  (j  <>  badbaud)  and  (j  <>  badbaud2)  then 
begin 

case  PHASEBITDIFF  of 
$01:  pbdl  :=pbdl+l; 

$02:  pbdl  :=pbdl+l; 

$04:  pbdl  :=pbdl+l; 

$08:  pbdl  :=pbdl+l; 

$03:  pbd2  : =pbd2+l; 

$05:  pbd2  :=pbd2+l; 

$06:  pbd2  :=pbd2+l; 

$09:  pbd2  : =pbd2+l; 

$0A:  pbd2  :=pbd2+l; 

$0C:  pbd2  :=pbd2+l; 

$07:  pbd3  :=pbd3+l; 

$0B:  pbd3  :=pbd3+l; 

$0D:  pbd3  :=pbd3+l; 

$0E:  pbd3  :=pbd3+l; 
end;  (*case  PHASEBITDIFF*) 


C*  now  count  the  total  number  of  phase  decoding  errors  *) 

TOTPHASERRS  :=  TOTPHASERRS  +  PHASEBITDIFF  and  $01; 
TOTPHASERRS  :=  TOTPHASERRS  +  (PHASEBITDIFF  and  $02) 
shr  1; 

TOTPHASERRS  :=  TOTPHASERRS  +  (PHASEBITDIFF  and  $04) 
shr  2; 

TOTPHASERRS  :=  TOTPHASERRS  +  (PHASEBITDIFF  and  $08) 

shr  3; 

end; 

(*  now  differentially  decode  the  magnitudes  of  the  tones  to  get  the 
fourth  and  last  bit  in  the  symbol  *) 

if  ( XNEWMAG  >  1.  5*X0LDMAG)  or  (XNEWMAG  <  2*X0LDMAG/3) 
then  XMAGBIT  :=  $01 
else  XMAGBIT  :=  $00; 

if  (RNEWMAG  >  1.  5*R0LDMAG)  or  (RNEWMAG  <  2*R0LDMAG/3) 
then  RMAGBIT  :=  $01 
else  RMAGBIT  :=  $00; 

if  (j  <>  badbaud)  and  (j  <>  badbaud2)  then 
begin 

T0TMAGERRS  :=  TOTMAGERRS  +  (XMAGBIT  xor  RMAGBIT); 
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numbits:  =numbits+4; 
end; 

(*  assign  colors  to  the  text  that  is  in  error  *) 

if  PHASEBITDIFF  >  0  then 

textcolor(138);  (*1. green  -  phase  error*) 

if  RMAGBIT  <>  XMAGBIT  then 

textcolor( 142);  (*yellow  -  mag  error*) 

if  (RMAGBIT  <>  XMAGBIT)  and  (PHASEBITDIFF  <>  0)  then 
textcolor( 140);  (*1.  red  -  dual  error*) 

if  (RMAGBIT  =  XMAGBIT)  and  (PHASEBITDIFF  =  0)  then 
textcolor( 15); 

(*  now  put  all  the  bits  together  and  color  the  errors  *) 

(*fill  TEMPBYTE  with  two  symbols*) 
if  f rac( SYMBOLCOUNT  /  2)  *  0. 5  then 
begin 

hue  :  =  textattr; 

XTEMPBYTE  :=  ((XPHASEBITS  or  XMAGBIT)  shl  4) 
RTEMPBYTE  :=  ((RPHASEBITS  or  RMAGBIT)  shl  4) 
end; (*  if  frac  *) 

if  (frac( SYMBOLCOUNT  /  2)  =  0.0)  then 
begin 

if  (hue  =  140)  or  (textattr  =  140)  then 
textcolor( 140); 

if  (hue  =  142)  and  (textattr  =  138)  then 
textcolor( 140); 

if  (hue  =  138)  and  (textattr  =  142)  then 
textcolor( 140); 

if  (hue  =  142)  and  (textattr  =15)  then 
textcolor( 142); 

if  (hue  =  138)  and  (textattr  =15)  then 
textcolor( 138); 

XTEMPBYTE  :=  (XPHASEBITS  or  XMAGBIT)  or 
XTEMPBYTE; 

XTEMPCHAR  :=  chr( XTEMPBYTE); 

RTEMPBYTE  :=  (RPHASEBITS  or  RMAGBIT)  or 
RTEMPBYTE* 

if  (RTEMPBYTE  =  $20)  and  (textattr  <>  15) 
then  RTEMPBYTE  :=  $5f; 

RTEMPCHAR  :=  chr( RTEMPBYTE); 

(*put  ascii  character  to  screen*) 

write( rtempchar); 

XTEMPBYTE: =0; 

RTEMPBYTE:  =0; 
end;  (*if  frac*) 

end;  (*for  decode  xmit  and  rec  data*) 

(*  now  calculate  the  means  and  variances  and  snrout  *) 
if  (j  <>  badbaud)  and  (j  <>  badbaud2)  then 
begin 

for  i:  =  1  to  8  do 
begin 

b:  =round( statmat[ i, 1] ); 
btot: =btot+b; 
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c:  =round(statoat[ i,2] ); 
ctot: =ctot+c; 
d:  ssro\ind(statmat[  i,3]  ); 
dtot:  =dtot+d; 

if  (i=l)  or  (i=4)  or  (i»5)  or  (i=8)  then 
begin 

for  count  :=  1  to  b  do 
begin 

big: =abs( recdata] ( 2*i)* 1 , count] ) 
+  big; 

ami:  *abs(recdata[ (2*i) .count] )+ 
ami; 


end; 

for  count  : «  1  to  c  do 
begin 

mbig:  *abs(recdata[  (2*i)-l+16 


count]  )+obig; 

msml:  «abs(recdata[ (2*i)+16, 
count] )+osml; 
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end; 

for  count  :=  1  to  d  do 
begin 

obig:  =abs(recdata[ (2*i)-l+32, 
count] )+obig; 

osml: =abs( recdata] ( 2*i)+32 , 
count] )+osml; 


end; 


end 

else 

begin 

for  count  :*  1  to  b  do 
begin 

big: =abs( recdatat (2*i) .count] 
+big; 

sral:=abs(recdata[(2*i)-l, count]  ) 
+sml; 

end; 

for  count  :  *  1  to  c  do 
begin 

obig:  =abs(recdata[ (2*i)+16, 
count] )+mbig; 

osml: =abs( recdata] (2*i)”l+l6, 
count] )+osml; 

end; 

for  count  :*  1  to  d  do 
begin 

obig:  *abs( recdata[ ( 2*i)+32 , 
count] )+obig; 

osml:  =abs(recdata[ (2*1) -1+32, 
count] )+osol; 

end; 

end; 

end; 

if  btot  >  1  then 
begin 

meanbig  :  *  big/btot; 
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meansml  :  ■  sml/btot; 
end; 

if  ctot  >  1  then 
begin 

mmeanbig  :  -  mbig/ctot; 
mmeansml  :  =  msml/ctot; 
end; 

if  dtot  >  1  then 
begin 

omeanbig  :  =  obig/dtot; 
omeansml  :  =  osml/dtot; 
end; 

for  i:=  1  to  8  do 
begin 

b: =round( s  t  atmat [ i , 1] ) ; 
c: =round(statmat[ i,2] ); 
d:  =round(statmat[  i,3] ); 
if  (i=l)  or  ( i=4)  or  (i=5)  or  (i=8)  then 
begin 

for  count  : =  1  to  b  do 
begin 

del:  =del+sqr(abs( recdata] (2*i)-l, count] )- 
meanbig)+sqr(abs( recdata] (2*i) .count] )-meansml); 
end; 

for  count  :=  1  to  c  do 
begin 

mdel:  =mdel+sqr(abs(recdata[ (2*i)- 1+16, count] )- 
mmeanbig)+sqr(abs(recdata[ (2*i)+16,co\uxt] ) -mmeansml); 
end; 

for  count  :  =  1  to  d  do 
begin 

odel:  =odel+sqr(abs(recdata[ (2*i)-l+32, count] )- 
omeanbig)+sqr( abs( recdata] ( 2*i)+32 , count] ) -omeansml) ; 
end; 

end 

else 

begin 

for  count  :  =  1  to  b  do 
begin 

del:  =del+sqr(abs( recdata] (2*i) .count] )- 
meanbig)+sqr(abs( recdata] (2*i)-l, count] )-meansml); 
end; 

for  count  :  =  1  to  c  do 
begin 

mdel: =mdel+sqr(abs( recdata] (2*i)+16, count] )- 
mmeanbig)+sqr( abs ( recdata] ( 2*i) -1+16 , co;mt] ) -mmeansml); 

end; 

for  count  :=  1  to  d  do 
begin 

odel: =odel+sqr(abs( recdata] (2*i)+32, count] )- 
omeanbig)+sqr( abs( recdata] ( 2*i) - 1+32 , count] ) -omeansml); 

end; 

end; 

end; 
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(*  use  a  complex  magnitude  in  statistical  calculations  only  if  it  had 
more  than  one  occurrence  *) 
if  (btot  >1)  then 
begin 

varx:  =de 1/ ( 2*btot ) ; 

snrbig:  *sqr(oeanbig) /varx; 

snrsml:  =sqr(oeansml)/varx; 

xmagbig:  =meanbig/cos(22. 5*pi/180. 0); 

xmagsml: =meansml/sin(22.  5*pi/180.  0); 

snravg:  =sqrC(xmagbig+xmagsml)/2)/varx; 

totsnr:  =totsnr+10*ln( snravg) /ln( 10. 0); 

bj:=bj+l; 

end; 

if  (ctot>l)  then 
begin 

mvarx: =mdel/(2*ctot); 
msnrbig:  =sqr(mmeanbig) /mvarx; 
msnrsml: =sqr(mmeansml)/mvarx; 
mxmagbig: =mmeanbig/cos(22. S*pi/180. 0); 
mxmagsml:  =mmeansml/sin( 22. 5*pi/180. 0); 
msnravg:  =sqr ( ( mxmagb ig+mxmagsm 1 ) / 2 ) /mvarx; 
mtotsnr:  =mtotsnr+10*ln(msnravg)/ln( 10.  0); 
cj: — cj+1; 
end; 

if  (dtot  >1)  then 
begin 

ovarx:  =odel/( 2*dtot); 
osnrbig:  =sqr(omeanbig) /ovarx; 
osnrsml:  =sqr( omeansml)/ovarx; 
oxmagbig:  =omeanbig/cos(22. 5*pi/180. 0); 
oxmagsml:  =omeansml/sin(22.  5*pi/180.  0); 
osnravg:  =sqr( (oxmagbig+oxmagsml)/2) /ovarx; 
ototsnr:  =ototsnr+10*ln(osnravg)/ln( 10. 0); 
dj: =dj+l; 
end; 
end; 

end;  (*for  )  :=  1  to  numbauds*) 
wr iteln( output ) ; 

writeln( output, 'The  overall  inner  SNR0UT  is  ' ,(totsnr/bj): 8: 3, 

’  db')j 

writeln( output, 'The  overall  middle  SNR0UT  is  ' ,(mtotsnr/cj):  8:  3, 
’  db'); 

writelnCoutput, 'The  overall  outer  SNR0UT  is  ' ,(ototsnr/dj):  8:  3, 
’  db’); 

wr iteln( output ) ; 

writeln( output, 'Total  phase  decoding  bit  errors  = 

T0TPHASERRS:  5:  0); 

writeln(output, '(' ,pbdl, '  symbols  with  one  bit  phase  decoding 
error)'); 

writelnCoutput, '(' ,pbd2, '  symbols  with  two  bits  phase  decoding 
error)'); 

writelnCoutput,  C  >pbd3,  symbols  with  three  bits  phase  decoding 
error)’); 


66 


end. 


writeln( output ,'Total  magnitude  decoding  bit  error*  « 
USHtort!’ out  of  ' .numbits,'  bits  transmitted' 


close(recdat); 
close(xmitdat); 
close( output); 


« 
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